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A collection of results is presented relevant for the analysis of azimuthal asymmetries in inclusive 
two-hadron production at BELLE. The aim of this overview is to provide theoretical ingredients 
necessary to extract the Collins effect fragmentation function. The latter arises within the CoUins- 
Soper factorization formalism, which describes both the transverse momentum and Q^ dependence 
of the cross section and its angular dependences at low and moderate transverse momentum. Since 
the Collins effect is not the only source of angular dependences, a discussion of various other effects 
is included. This concerns higher twist contributions, photon-^-boson interference effects, radiative 
corrections, beam polarization and weak decays. Furthermore, different frames, transverse momen- 
tum weighting and ratios of asymmetries are discussed. These issues are all of relevance for the 
5^ , unambiguous measurement of the Collins effect. 
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t"-- , I. INTRODUCTION 



(~| . The BELLE experiment at KEK in Japan measures with very high luminosity the process of electron-positron 

Q^' annihilation in collisions of 8.0 GeV electrons and 3.5 GeV positrons. The center of mass energy is selected to be 
on-resonance of the T(4S') meson, which has a mass of 10.5800 ± 0.0035 GeV and which decays more than 96% of 
(Vj ' the time into BB meson pairs. The main aim of the measurements of B and B decays is to study CP violation. 
r^ \ Besides this goal, there are other interesting studies that can be performed at BELLE^ and for which also the 
off- resonance data are useful (which in the case of BELLE are taken 60 MeV below resonance). The acquisition 
of off-resonance data is mostly used for background studies, but is also of interest for physics studies that are not 
^ . 6-quark specific. This overview discusses such a case, namely the study of azimuthal asymmetries in the inclusive 
00 ' production of two almost back-to-back hadrons, e+e" -^ hi h2X. There are several effects that can cause such 
f^ , asymmetries and one would like to disentangle them in order to isolate perturbative from nonperturbative effects. 
■^ ' In the latter type of effects quark spin is expected to play a nontrivial role via the so-called Collins effect [l| . By 
^^ \ measuring azimuthal asymmetries it may be possible to obtain trustworthy quantitative knowledge on this type 
Tij- ■ of effect, for instance from the BELLE data. Recently, the first results from BELLE were published [2'| and the 
CD \ purpose of this overview is to discuss the theoretical aspects of such type of study. 

OO i The Collins effect was first discussed in the context of semi-inclusive deep inelastic scattering (SIDIS) of leptons 

^7^ ' off transversely polarized protons, as a means to access transversity \^. Transversity [3| describes the extent to 

which quarks are transversely polarized inside a proton that is polarized transversely to the probing particle, which 

in the case of SIDIS is a virtual photon. The Collins effect describes the angular asymmetry in the distribution 

^ . , of hadrons produced from a transversely polarized fragmenting quark. Via this effect the transverse polarization 

j^ ' of the struck proton results in an asymmetric distribution of final state hadrons. If sufficiently large, the Collins 

effect would thus allow for a measurement of transversity and subsequently of the tensor charge, the fundamental 

charge that can only be measured through transversity. 

The first nonzero Collins effect asymmetry in polarized SIDIS has been observed by the HERMES experiment 
[J] (using a deuteron target the COMPASS experiment obtained a result consistent with zero ^], presumably due 
to cancellations between proton and neutron contributions). The HERMES result indicates that both transversity 
and the Collins effect are nonzero. For an extraction of transversity from those SIDIS data a separate measurement 
of the Collins effect fragmentation function needs to be performed. This can only be done in the process e^e~ — *■ 
hi h2 X [6| and motivates the BELLE efforts concerning the measurement of this process. Some earlier attempt 
to use LEPl data has been undertaken ,7], but without study of systematic effects and still remains preliminary. 
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^ We have chosen to focus on BELLE, but the observables to be discussed can of course also be studied at other e+e~ colliders. 



Moreover, as will be discussed, it is likely that the Collins effect asymmetry in e+e~ ^ hih2X has a powerlike 
fall-off behavior with energy, which would favor an extraction at BELLE over LEPL 

We will study various effects that could lead to azimuthal asymmetries in the process of interest, e^e'^ — > 
hi h2X, which besides the Collins effect, include electroweak j-Z interference effects, beam polarization effects 
and radiative corrections. In order to arrive at an unambiguous interpretation of the data the magnitude and 
scale dependence of the various effects need to be estimated, of course to the extent to which that is possible from 
first principles. The effects will not be treated simultaneously; combinations of effects will only be considered 
when the analysis requires it. As a rule we will ignore effects that are expected to be smaller than a permille, such 
as Z-Z contributions or beam polarization in combination with j-Z interference. From the possible contributions 
considered the higher-twist effects are the least known, but unfortunately not necessarily below the percent level 
at BELLE. In cases where no reliable estimate can be given, such as for twist-4 effects (©(A^/Q^)), additional 
observables or checks may need to be considered to further exclude competing interpretations of the Collins effect 
asymmetry. But despite some uncertainties the isolation of the Collins effect contribution does seem feasible given 
the possibilities that BELLE offers. The purpose of this overview, which contains several new aspects, is to assist 
and facilitate this endeavor. The realistic prospect of extracting the Collins effect fragmentation function and with 
it transversity, of which a first result has recently been obtained [8| , makes the study of azimuthal asymmetries 
at BELLE well worth the effort. 

This overview consists of the following sections and subsections: 

im General angular dependence 

mil Two-particle inclusive cross section 

nil Al Frames 

IIIIBI Structure functions 

IIVI Leading order cross section 

IIV Al Integration over transverse photon momentum 

IIVBI Unintegrated cross section 

IIV CI Weighted cross sections 

IIVDI Estimate of the Q|n-weighted Collins effect asymmetry 

IVl Universality of the Collins effect 

rvn Higher twist 

IVIII Electroweak interference effects 

IVIIII Jet frame asymmetry 

IIXI Scale dependence of the Collins effect asymmetry 

IIX Al CoUins-Soper factorization 

IIXBI Numerical study of the Q^ dependence of the Collins effect asymmetry 

IIX CI Comparison to tree level 

IIXDI Nonpcrturbative Sudakov factor from BELLE 

1x1 Radiative corrections 

IXII Weighted asymmetry beyond tree level 

Kill Ratios of asymmetries 

IXIIII Beam polarization 

IXIVI Weak decays background 

IXVI Summary 

II. GENERAL ANGULAR DEPENDENCE 

We consider e^e+ —>■ hi ft,2 A, where the two leptons (with momentum / for the e^ and V for the e'^) annihilate 
into a photon (or Z boson) with momentum q = I + 1' . This photon momentum sets the scale Q, where Q^ ^ q^, 
which is much larger than characteristic hadronic scales. Denoting the momentum of outgoing hadrons by Ph {h 
= 1, 2) we use invariants Zh = IPh ■ qjQ^ ■ We will consider the case of unpolarized leptons and hadrons, although 
in Sec. IXIIII we will turn to the issue of transverse beam polarization due to the Sokolov-Ternov effect. We will 
work in the limit where Q^ and P^ ■ q are large, keeping the ratios Zh finite. We will consider the case where the 
two hadronic momenta Pi and P2 do not belong to the same jet (i.e.. Pi ■ P2 is of order Q^). 

There is a considerable literature on angular correlations for three-jet events in electron-positron annihilation, 
e.g. [3, [10, [lH [12, [l3, H, m, H, [l3, [13] , but here we are Hmiting the discussion to two-jet events exclusively. 



dN / da \ da 



dfl \dzidz2d'^qrp J dzidz2dQd'^ Qrp 



In general, the differential cross section for the process e+e ^ hih2X can be written as (see for instance |14| ) 

Fi{l + cos^ 61) + F2(l - 3 cos^ 9) + Fg cos 
+F4 sin 26 cos (/) + F5 sin^ 6 cos 20 + Fg sin 6* cos 

+F7 sin 26* sin + ^8 sin^ 8 sin 20 + Fg sin 6* sin 0. (1) 

The functions Fi depend on the invariants Zh — 2Ph ■ q/Q^ and on q^ = Q^, the squared transverse momentum 
Qrp of the photon with respect to the two hadrons. The angles (j) and 9 are given in the lepton-pair center of mass 
frame or equivalently the photon center of mass frame. Precise definitions and explanations will be given below. 

If at high Q^ and Q^ collinear factorization of the cross section is considered, then at tree level (zeroth order in 
as) only Fi, F3 will receive nonzero contributions (F3 only from 7-Z interference), at first order in as Fi, . . . , Fg 
receive contributions and at second order all Fi are nonzero. All this is assuming no transverse beam polarization 
is present. As said, the complication of transverse beam polarization will be considered in Sec. IXIIII 

Below we will study the differential cross section da / dzidz2dQ.d^ qrp in much detail, also at lower values of Q'^ 
where collinear factorization is not the appropriate framework. To set the notation we will first look at the cross 
section expression in terms of the hadron tensor. 

III. TWO-PARTICLE INCLUSIVE CROSS SECTION 

The square of the amplitude for e^e+ -^ hi h2 X can be split into a purely leptonic and a purely hadronic 
part, 

|Mp = ^F^.F''^ (2) 

with the helicity-conserving lepton tensor (neglecting the lepton masses) given by 

Lfj,^{l, I') = 2lfj_ll + 21J'^ - Q'^g^u- (3) 

For the case of two observed hadrons in the final state, the product of hadronic current matrix elements is written 
as 

H^.,{Px;Pi;P2) = {0\MO)\Px;Pi;P2){Px;Pi;P2\Jumo). (4) 

The cross section for two-particle inclusive e'^e" annihilation is given by (including a factor 1/2 from averaging 



over initial state polarizations) 



L^.W-', (5) 



d^Pid^P2 4g6 "" 

with 

W^.(9;Pi;P2) = J (2^)3^po S\q - Px - Pi - P2)i?p.(Px;Pi;P2). (6) 

A. Frames 

For the calculation of the hadron tensor it will be convenient to define lightlike directions using the hadronic 
(or jet) momenta. The two hadronic momenta Pi and P2 can be parameterized using the dimensionless lightlike 
vectors n+ and n_ (satisfying n^ ■ n_ = 1), 

P^^^n^_ + ^!l^nl, (7) 

V2 C1QV2 + 

P^ = —11^ n^ + ^ n!f , (8) 



M2 fj. , C2Q ^ 
C2Q\/2 72 + 

g^^-|n'i + -|<+g^, (9) 



Q^ with q^ = —Q'^. When Q"^ <^ Q^ one has Q = Q, Ci — zi and (2 — ^2 up to Q^/Q^ 



where Q^ — Q 
corrections. 

Vectors transverse to 71+ and ri_ one obtains using the tensors 



9t 



flV 
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(10) 
(11) 



The experimental analysis of the azimuthal asymmetries will usually not be performed in the frame in which 
the two hadrons are coUinear, for which the above (Sudakov) decomposition into lightlike vectors and transverse 
parts is most suited. Instead it is much more common to consider angles in the lepton-pair center of mass frame 
or equivalently, the photon rest frame. In this case there is still freedom to select which momentum (or linear 
combination of momenta) is used to define the z axis, or equivalently, what determines the perpendicular plane 
in which the azimuthal angles lie. A few choices are common, such as the Gottfried-Jackson frame and the 
Collins-Soper frame. 

For the most part of this overview the frame to be employed will be the e+e^-annihilation analogue of the 
so-called Gottfried- Jackson frame [T^]. We will first give the details of this frame. Later on we will also consider 
the analogue of the Collins-Soper frame [20| and a frame where the jet or thrust axis is used to fix the basis 
(explained separately in Sec. IVIII[) . We will alternate between these different basis sets depending on what is 
most convenient for the analysis. 

In order to expand the hadron tensor in terms of independent Lorentz structures (parameterized by structure 
functions), it is convenient to work with vectors orthogonal to q. A normalized timelike vector is defined by q and 
a normalized spacelike vector is defined by P^ = P/* — {Pi ■ q/q'^) q^ for one of the outgoing momenta, say P2, 






P2-q ^ 



P^ (7*' 

Z2Q Q 



(12) 
(13) 



This choice of frame is the analogue of the Gottfried-Jackson frame often employed for the Drell-Yan process. 
This means that in the lepton-pair center of mass frame, hadron 2 is moving along the z direction (see Fig. [T]). 
In general, in this frame hadron 1 will have momentum components orthogonal to z and t. 




lepton plane (cm) 



FIG. 1: Kinematics of the anniliilation process in the lepton center of mass frame (the analogue of the Gottfried-Jackson 
frame) for a back-to-back jet situation. P2 is the momentum of a hadron in one jet. Pi is the momentum of a hadron 
belonging to the other jet. 



Vectors orthogonal to z and i are obtained with help of the tensors 

M-' + z^'z", 
1 

e ^2p<la 



gT^g'" 






1 



{P2 ■ q) 



(14) 
(15) 



Since we have chosen hadron 2 to define the longitudinal direction, the momentum Pi of hadron 1 can be used 
to express the directions orthogonal to t and z. One obtains P^j_ — g^ Pi^, . We define the normalized vector h^ 
= P^j_/\Pi±\ and the second orthogonal direction is given by €^^h^. 

Note that the transverse tensors in Eqs. (jlOp and (|lip are not identical to the perpendicular ones defined above 
if the transverse momentum of the outgoing hadron 1 does not vanish. The lightlike directions can easily be 
expressed in t, z and h, namely, up to corrections of order Q\IQ^ , 



n_ 



1 

1 

71 



[i> + z^], (16) 

(17) 



t> - z^ + 2 % h^ 



This shows that the differences between g^^ and gl^ are of order 1/Q. Especially for the treatment of azimuthal 
asymmetries subleading in l/Q fSec. IVip . it is important to keep track of these differences. We will see however 
that taking transverse momentum into account does not automatically lead to suppression. 

Thus far we considered two sets of basis vectors, the first set constructed from the two hadron momenta (Pi and 
P2), the second set from the photon momentum (g) and one of the hadron momenta {P'l). The respective frames 
where the momenta P\ and P2, or g and P2, are coUinear are the natural ones connected to these two sets. One 
can go from one frame to the other via a Lorentz transformation that leaves the minus components unchanged 
[2l|. In the first frame q has a transverse component g^^, in the second P\ has a perpendicular component Pi^. 
We will therefore sometimes refer to them as the "transverse basis" and the "perpendicular basis" , respectively. 
Up to corrections of order Q\IQ^ qij^ and P^j_ are related as follows: 

P/Y = -zig^ = ziQT/j^. (18) 

Azimuthal angles will lie inside the plane orthogonal to t and z. In particular, cj)^ gives the orientation of 
Z_L, where l'^ denotes the normalized perpendicular part of the lepton momentum l^^. The angle 4>i is between 
h ex Pi± and l±. More specifically, 

f_L • a_L = -|a_L|cos0a, (19) 

f^T^^f^"'^!^ = \a±\sin(f>a, (20) 

for a generic vector a. The convention for the epsilon tensor used is e'^^^'^ = 1. 

Sometimes it may be convenient to choose a different (rotated) set of basis vectors in the lepton center of mass 
frame, the Colhns-Soper frame. For comparison let us denote the basis vectors of the Gottfried- Jackson (GJ) 
frame for e+e^ —f hih2X as: 



-GJ 



"-GJ 



Q 5M 



P2-q 
2Q_ 

sQt 



P2^ (22) 

(Pi • ZGj)P^ - (P2 • ZGj)P^) , (23) 



where s = (Pi +P2)^. Here xgj corresponds to h and yoj is defined implicitly by requiring a right-handed basis. 
The basis for the CoUins-Soper (GS) frame for e+e^ -^ hih2X is defined as |22 |: 



gM 



^((Pi.<z)P2^-(P2-g)Pf), (25) 



"■GS 



^GS 



SQ 

2Q 



sQtQ 



^({P,-q)P!^ + {P2-q)P^). (26) 



Note that we are keeping terms of order Q\IQ^., but not order M^ jQ"^ . Throughout this paper we will neglect 
target mass corrections. 




lepton plane (cm) 

FIG. 2: Kinematics of the annihilation process in the lepton center of mass frame, the analogue of the Collins-Soper frame. 

In the lepton-pair center of mass frame the z axis now points in the direction that bisects the three- vectors P2 
and —P\ (see Fig. [5|). One finds that in the limit Qt ^^ 0: zcs a-nd zqj coincide and also xcs and xqj. When 
Qt 7^ they differ only by a rotation: 



zgj = COS/3 Zcs + sin/3 xcs, 
Xqj = -sin/3 zcs + cosP xcs, 



(27) 
(28) 



where 



cos /3 = — , sm /3 = ^^ . 
Q Q 



(29) 



We will also refer to the Collins-Soper frame vectors i, z,x,y as a perpendicular basis. 

In the cross sections one will encounter the following functions of ?; = P2 ■ 1/ P2 • Q ~ ^ Z? i which in the lepton- 
pair center of mass frame equals ?; = (1 -f cos9*)/2, where 6* is the angle of z with respect to the momentum of 
the incoming lepton I: 



Aiy) 



1 



B{y) = y{i-y) =' ^ sin^ 



C{y) - (l-2j;) 



D{y) - Vvi^-y) 



cosr, 

isin6 



(30) 

(31) 
(32) 
(33) 



In the Gottfried- Jackson frame (Fig.[T]) 9* is called 82 and in the ColHns-Soper frame 9 (Fig. [2]). 

The cross sections are obtained from the hadron tensor after contraction with the lepton tensor. The lepton 
tensor for unpolarized leptons expressed in the lepton center of mass is given by {l'^ = l'^/{Qy/yiT^^y))): 



LI"" 



Q' 



-2A{y)gY + AB{y)z>^z'^ - AB{y) (l^Jl + \ g^) - 2C{y)D{y) 2^%^ 



(34) 



For later use, the contractions of specific tensor structures in the hadron tensor are given in Table HI 
In the e+e" center of mass frame dP PidP P2 / Pi P2 = {dzi/zi){z2Q^dz2/A)cPPi±dfl2, such that 



da'^'^*^ 



dzidz2dD,d'^qrp 16(5^ 
where dfl — 2dyd(j)^, with <j)^ giving the orientation oi l±. 



ziZ2L^,W^'', 



(35) 



TABLE I: Contractions of the lepton tensor Lp^ with tensor structures appearing in the hadron tensor. 



iM.™''7(4Q' 



-9T [h-y + y') 

aY'bY - (a± • fax) ffx" ^y (1 ^ v) I'^J-I l^'J-l cos(0a + 0b) 

£("«"* -(l-2y)^y(l-j;)|ax| cos0. 



B. Structure functions 

The hadron tensor W^"^ can be expanded in terms of independent Lorentz structures which leads to a parame- 
terization in terms of structure functions Wi. Ignoring lepton polarization and 7-Z interference, the most general 
decomposition consists of four structure functions. Due to the similarity of the process e^e~ —^ hi h2 X with the 
Drell-Yan process, we will emplo y si milar notation here as used for the latter process, i.e. we follow the notation 
of Lam & Tung [II,!!!, Collins [H], and Argyres & Lam ^. 

The structure functions in the leptonic center of mass frame (or rather the perpendicular basis) are defined as 

W""" = -g^^Wr + zf^z^WL - z^'^x^'^Wa - (x^'^x''^ - x2c,'['')PFaa, (36) 

such that W^^ = — (2Wt + Wl)- The structure functions Wt. l. A, aa are associated with specific polarizations 
of the photon [23;]: Wt = W^^'^, Wl = W°'°, Wa = (M^"'^ + W'^'")/V2, and Waa = W^'-^ where the first and 
second superscripts denote the photon helicity in the amplitude and its complex conjugate, respectively. In terms 
of these structure functions the cross section becomes 



dzidz2dn,(Pqrp AQ"^ 

or 

dN f da \~^ da 3 1 



(37) 



[Wt (1 + cos- 



2, 



dO. \dzidz2d^qrp J dzidz2dVt*d'^qj, Stt 2Wt + Wl 

+ M^l(1 - cos^ r) + Wa sin2r cos^* + Waa sin^ 0* cos 20*] . (38) 

Here d*,<j)* indicate the polar and azimuthal angle in the e+e^ center of mass frame (which for the BELLE 
experiment is not the lab frame), such that in the Gottfried- Jackson frame 9* = O2 and 0* — 4>i, the angles we 
have used before (see Fig. [T|) . 

Another standard notation for the angular dependences in the Drell-Yan process can be employed here as well: 



dN 3 1 



1 + A cos"' 9* + fi sin 20* cos </>* -t- - sin^ 0* cos 20* 



dQ 47r A + 3 

Expressing these parameters in terms of the structure functions W one has: 



(39) 



Wt- Wl Wa 2Waa . . 

Wt+Wl' ^ Wt + Wl' "" Wt + Wl' ^ ' 

One can transform from the CS frame to the GJ frame by using the following transformation matrix: 

1_1.2 _3 32 \ /;^ 




1 
GJ \ ^ -f- ^ ' 2'- / \ " / CS 



(41) 



V 



where p = Qt/Q and 



A = 1 + p2 + -p^X + p^ _ _p2j,^ 



(42) 



The transformation from the GJ frame to the CS frame is the same, but with the replacement p -^ —p. 

In the discussion of the contributions to the various structure functions given below, the structure function 
VFaa and the analyzing power v of the cos 2(j) asymmetry will receive particular attention due to its Collins effect 
contributions. 



IV. LEADING ORDER CROSS SECTION 

In this section we investigate the 4> dependence that arises in leading order in as and 1/Q in the cross section of 
the process e+e" ^ hi h2 X differential in the transverse momentum Pi± = —ziqj,. The cross section involves 
products of frag mentation functions, which unlike the ordinary collinear functions include transverse momentum 
dependence [27| . The idea that such "intrinsic" transverse momentum will give rise to power suppression turns 
out not to be true, even though this was true in the pioneering studies |28l . |29| on azimuthal dependences due to 
intrinsic transverse momentum. Nontrivial quark spin effects, which require nonzero partonic transverse momenta, 
can arise at leading power. One such effect is the Collins effect, which gives rise to a cos 20 asymmetry. This was 
first pointed out in Refs. [6|, |30| and in this section we will repeat the essentials. 

A. Integration over transverse photon momentum 

Although we are interested in the two-hadron inclusive cross section differential in q^p , we will first consider the 
case of integration over the transverse momentum of the photon. At tree level one needs to calculate the diagram 
shown in Fig. [31 It depicts the squared amplitude of the process in which the photon produces a quark and 




FIG. 3: Factorized diagram contributing to e e annihilation in leading order. There is a similar diagram with reversed 
fermion flow. 

an antiquark, which subsequently fragment independentl y in to the hadrons hi and /12, respectively. The quark 
fragmentation correlation function A (Pi; k) is defined as [27| : 



A,,{Pi;k) = J2 f 
X •' 



(2^ 



e*-(0|V'.(x)|Pi;X)(Pi;X|V',(0)|0), 



(43) 



where k is the quark momentum and an averaging over color indices is left implicit. It is also understood that 
appropriate path-ordered exponentials should be included in order to obtain a color gauge invariant quantity, 
cf. e.g. [3l|. In Sec. |V] we will briefly address the universality issue that arises from the proper gauge invariant 
definition. 



The above matrix element as function of invariants is assumed to vanish sufficiently fast above a characteristic 
hadronic scale {0{M)) which is much smaller than Q. This means that in the above matrix elements fc^, A: • Pi <C 
Q^. Hence, we make the following Sudakov decomposition for the quark momentum k: 

fcEE-!-^n_ + -i —^n+ + kr~-Pi + kr. 44) 

zV2 Z1QV2 ^ ^ z ^ 

The Dirac structure of the quark correlation function can be expanded in a number of amplitudes, i.e., func- 
tions of invariants built up from the quark and hadron momenta, constrained by hermiticity and parity. In the 
calculation of the cross section integrated over the transverse momentum of the photon, at leading order we only 
encounter the integrated correlation function / dk^ (Pkx^iPi; k), which is a function of k^ only. At leading twist 
this leaves only one possible Dirac structure: 



A(z) EE J f dk+(fkT A(Pi;fc) 



k-=p-/z "^Pl 



-Di{z)ri- (45) 



The function -Di(z) is the ordinary unpolarized fragmentation function. 

For the fragmentation of an antiquark most things are analogous to the quark fragmentation. The main 
difference in the present case is that the role of the + and the — direction is reversed. We will denote the 
antiquark correlation function by A{P2;p): 

A,,{P,;p) = E/ (^ ^"''" m,mP2;X){P2;Xm^)\0). (46) 

Similarly, 



A{z)= ^ j dp- d^p^ A{P2;p) 



p+=p:^/z *^2 



4Po 



-Di{z)^2- (47) 



The four-momentum conservation delta-function at the photon vertex is written as (neglecting 1/Q^ contribu- 
tions) 

5\q -k-p) = 5{q+ - p+) S{q- - k') S^p^ + k^ - q^), (48) 

fixing P2 /z = P^ = q^ = P2 / '^'2 and Pi j z = k^ = q^ = P^ /zi. Eq. (|48p shows why only the fc+ and 
^"-integrated correlation functions are relevant. 
The hadron tensor as function of q^^ is given by 

W^" = 3 / dp-dk+d^pj,d^kj, S^{pt + fcy - Qt) Tr (A{p)j'' A{k)j'') . (49) 

The factor 3 originates from the color summation. We have omitted the flavor indices and summation; furthermore, 
there is a contribution from diagrams with reversed fermion flow, which results from the above expression by 
replacing fi -i^ v and q — > —q (and in the end in a summation over flavors and antiflavors). Note that the quark 
and antiquark transverse momentum integrations are linked, unless one integrates over qrp. 

After integration over the transverse momentum of the photon (or equivalently over the perpendicular mo- 
mentum of hadron 1: Pi±^ = —ziqj,), the integrations over krp and Pj, in the hadron tensor in Eq. (j49p can be 
performed resulting in^ 

d\^ W^"^ = _il^e^gr^i^i: (50) 

We have now included the summation over flavor indices and e^ is the quark charge in units of e. The frag- 
mentation functions are flavor dependent and only depend on the longitudinal momentum fractions, i.e. -Di-Di = 

D1[zi)Dl{z2). 



^ Eq. II50I I corrects Eq. (75) of Ref. jj| by an additional factor (Z1Z2) 
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From this hadron tensor one arrives at the foUowing expression for the cross section at leading order in as and 

i/Q 



da{e+e ^ hih2X) 3a^ ,-^ 

- = -TrT My) } e^DiDi 



dzidz2d£l 



a. a 



(51) 



B. Unintegrated cross section 



Now we turn to the cross section differential in the transverse momentum. In this case the correlation function 
A only integrated over fc+ needs to be considered. It can be parameterized in terms of transverse momentum 
dependent (TMD) fragmentation functions [32, 123] 



A(z,M = 



1 

4^ 



dk+ A(Pi;fc) 



M- 



k-=P~ /z, fci 



^^{DMkD 



4R 



Ml 



H, 



-(z,4)''^'^^^^i' 



M( 



(52) 



where we only display the fragmentation functions that are relevant for unpolarized hadron production. The 
functions Di and H^ yield contributions to the cross section that are of leading order in 1/Q. After integration 
over fey the term with H^ drops out and the first term reduces^ to the expression in Eq. (|45l) . Strictly speaking, 
the TMD fragmentation functions depend on z (the lightcone momentum fraction z = P^ /k^ of the produced 



21.2 



Here k'rp = 



-zkj' is the transverse momentum 



hadron with respect to the fragmenting quark) and on fc^^ — z^k 

of the hadron in a frame where the quark has no transverse momentum. In order to switch from quark to hadron 

transverse momentum a Lorentz transformation leaving k~ and P{~ unchanged needs to be performed 21, 27]. 

The so-called Collins effect function H^ implies a correlation between the transverse polarization direction of 
the quark and the transverse momentum direction of the unpolarized hadron it fragments into [l| . The Collins 
effect correlates the azimuthal angle of the transverse spin of a fragmenting quark with that of the transverse 
momentum of the produced hadron (both taken around the quark momentum) , via a sin (j) distribution of their 
difference angle 4>. Therefore, the distribution of final state particles contains information about the spin direction 
of fragmenting quarks. In this sense it is the strong interaction analogue of the self- analyzing property of weak 
decays. 

The presence of the Dirac matrix cr^i^ shows that the Collins effect is a chiral-odd (a quark-chirality flip) state; 
an interference term between opposite chirality states of the fragmenting quark. The function H^ is also often 
referred to as 'time-reversal odd' fragmentation function, due to its behavior under time reversal. It does not 
imply a violation of time reversal symmetry though. For a detailed discussion cf. Ref. 31]. 

The cos 2(f) asymmetry to be discussed below depends on a product of two Collins effect fragmentation functions 
H^. It is an azimuthal spin asymmetry in the sense that the asymmetry arises due to the correlation of the 
transverse spin states of the quark-antiquark pair. On average the quark and antiquark will not be transversely 
polarized, but for each particular event the spins can have a transverse component and these components will 
be correlated via the photon polarization state, which in turn is determined by the lepton direction. Due to the 
Collins effect the directions of the produced hadrons are correlated to the quark and antiquark spin and hence, 
to the lepton direction. This correlation does not average out after summing over all quark polarization states. 

As a transverse spin state is a helicity-flip state, one deduces that the asymmetry arises from the interference 
between the photon helicity ±1 states (along the quark-antiquark axis) and hence contributes to Waa and to i'. 
Such a helicity-flip contribution can also arise from quark mass terms, but those are power suppressed and do not 
lead to an azimuthal dependence (the 9 dependence will be identical though). 

The Collins effect was the main reason for studying azimuthal asymmetries in the BELLE data, since it shows 
up at leading order (in both as and 1/Q) in an azimuthal cos2(^ asymmetry in the differential cross section for 
unpolarized e'^ e^ —^ hi h2 X [a|: 



do-(e^ 



dzidz2dfld-^qj 



hih2X) _ 3c^ 
" Q^ 



n^2 



A{y) T {DiDi\ 



^ The relation between TMD and ordinary collinear fragmentation functions is not trivial beyond leading order in o.s 
on the analogous problem for distribution functions we refer to Ref. [34l . |35(I . 



For a discussion 
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+B{y) cos 201 T 
where we use the convohition notation 



\ Zi fX' K^ lb' p^ — K^ ' Py' 



Hint 



T[DD] =^el f (fkrp (fpr^. 6^{prp + krp - qj.)D''{zi, zlk^)D''{z2, zIp^). 



(53) 



(54) 



The angle (f>i is the azimuthal angle oi h — x, see Fig. [T] So we find that in the lepton-pair center of mass frame: 

-1 



dN 

'dn 



da 



da 



dzidz2 d^ Qj, J dzi dz2 dfld^Qj 



—- { (1 + cosH2)T [DiDi] 



+ sin^ 02 cos 201 J^ 



(2h-krp h 



' P'T' K'Y' ' P'T' 



Ml Ma 



T[DiDi] 



This shows that 



Wl ^ Wa ^ 0, 



Waa 



T 



\ A it' rZi-p ft' P'T' ' fC T"! ' P'T' 



M1M2 



(55) 



(56) 

(57) 

(58) 



or equivalently, that at tree level A = 1, /i = and 



T 



v = 2- 



(2h-krp hpj, — krpPj, j H^H^ 



MiM2T\DiDi 



(59) 



We emphasize that measuring v does not involve a measurement of the polarization of the produced hadrons nor 
of the incoming leptons. Also, the result is not suppressed by factors of 1/Q, in contrast to the cos 20 asymmetry 
discussed by Berger [2^, which is l/Q^-suppressed. Eq. ([55)1 applies to the GJ frame, but upon neglecting Q\IQ'^ 
power suppressed terms it is the same expression in the CS frame, as can be seen using Eq. (|^T]) . 

We will often assume Gaussian fe^-dependence of the various functions, since in that case the convolutions can 
be explicitly evaluated. Eq. (|53p then becomes 



daie^e -^ hih2X) 3a' 



dzidz2dQd'^qrj 



= -^ Q{Qt\ ^) E ^a ^(y) DI{Z,)DI{Z2) 



+B{y) cos 201 



QrpR 



MiNhRiRi 



Ht^{z^)Ht%Z2 



where R^ = RlRl/iRl + Rl) and 



Diizi^k'^') = i?i(zi)i?f exp(-i?ffc^)/7rzf = i?i(zi)^(|fcj,|; i?i)/z^. 



(60) 



(61) 



and similarly for Di, H^, H^ with obvious replacements (for details cf. Ref. [33]). For simplicity we have assumed 
the same Gaussian width for Di and H^, which is not expected to be realistic. Later on we will drop this 
assumption. Rather, we will often assume Ri = R2 = R which means equal widths for H^ and Hj^, Riu = 
R2U — Ru which means equal widths for Di and Di, and moreover. Mi = A/2 — M, which should be a reasonable 
assumption when the two produced hadrons are charged pions. This leads to 



da{e^e 



dzidz2dVLd'^q 



hih2X) _ 3a2 ^ 



A{y) g{QY-Rj2) Dl{zi)Dl{z2) 



B{y) cos 201 ^ GiQT-.Rm Ht^{z^)Ht^{z2) 



Q\ 



(62) 
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C. Weighted cross sections 

The expressions in the previous subsection contain convolutions, which are not the objects of interest, rather one 
wants to learn about the fragmentation functions depending on z and fcj^. This may not be possible without further 
assumptions about the type of fc-r. dependence, such as assuming Gaussian transverse momentum dependence. As 
a way out, it has been suggested [6|, |36[ to consider specific integrated, weighted asymmetries that probe instead of 
the full transverse momentum dependence, the fej^-moments of the functions. These so-called transverse moments 
are defined as: 

F(")(zi) =|d^fc^ [^^^ nzu^T\ (63) 

for a generic fragmentation function F . In particular, the first transverse moment of the Collins fragmentation 
function 

H^^'^ (z) = z'J d'kr ^ Ht {z, z- k%) (64) 

has been considered frequently in the literature. 

In Sect. lIV"S1 we have presented the hadron tensor and cross section integrated over transverse momentum of the 
photon. A number of structures averaged out to zero, which are retained when the integration is weighted with an 
appropriate number of factors of q^. By constructing such weighted cross sections at tree level the convolutions 
become simply products of such fcj^-moments. The fej^-moments can be used in other processes where they also 
occur. 

To shorten the notation we define weighted cross sections as follows 

/MA [ j2 w d(T{e+e- -^ hih2X) 

i^)=JdlTW ^^^^^^^^^.^^ , (65) 

where W = W{Qt, <Pi)- In this way one finds: 

(1) = ^A{y)J2<^lDU^i)DU^2), (66) 

a,d 

(Qt) = ^2A{y)Y^el[M^D^^^''{z,)DU^,) + M,'DUzM^''iz2)), (67) 



/ Qt oj, \ 3a 



ViM^M^^^^^'^V ^ 'Q 






a,a 



COS201 = ^^B{y) y^et Ht'"" {z^) H ,' ' (z^). (68) 



The fcji-moment H^ that arises in the above e+e~-annihilation expression also appears in the Q^-weighted 
sin(0^ + (j)s) asymmetry in semi- inclusive lepton-hadron scattering, in that case multiplied by the transversity 
distribution function [37i | . This illustrates the purpose of considering such weighted cross sections. 

Next we will discuss an estimate of the tree level weighted expressions of interest for the Collins effect. In Sec. 
IXII we will address the effect of radiative corrections on such weighted asymmetries. 

D. Estimate of the Q|— weighted Collins effect asymmetry 

A natural question to ask is what one expects for the magnitude of the Collins effect asymmetry. This runs 
immediately into the problem that the Collins fragmentation function is a nonperturbative quantity that is at least 
as hard to calculate from first principles as the ordinary unpolarized fragmentation function Di . In this subsection 
we will discuss a rough estimate of the weighted cross section defined in Eq. (|68p . appropriately normalized. 

For an order of magnitude estimate of the weighted asymmetry, we consider the situation of the produced 
hadrons being a tt"*" and a tt^ and only consider up and down quarks. Furthermore, we assume D^^'" (z) — 
Df^^ (z), Df^'^ (z) = Di^'" (z) and neglect unfavored fragmentation functions like D'^^'^ (z), etc; and 
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similarly for the time-reversal odd functions H^,H^. These equalities seem quite safe for the Di functions on 
grounds of isospin and charge conjugation. The same assumptions might be non-trivial for the H^ functions and 
remain to be tested. Studies of the HERMES data seem to indicate that the unfavored Collins functions can be 
of the same magnitude as the favored ones [J] ■ It remains to be seen whether this also holds true at higher energy 
scales, when the multiplicity of hadrons in the final state is significantly higher and momentum conservation would 
not correlate the hadron momenta as much as at low multiplicities. 
With these assumptions we obtain 

\4M^ / Di{zi) Di{z2) 

where 

p. X B{y) cm SVC? 02 , , 

^(^) ^Ad)- l + cos^92- (^^^ 

An upper bound on the ratio H^ ^ '{zi)/Di{zi) could be obtained using the Soffer type of bound [3a.l3g| 

Ifc^l \Htiz, |fcj,|)| < Ah D,{z, Ifc^l), (71) 



which should hold for all \kj.\. It implies: 



' ^ ^ ^' < {\k^\){z)/{2M,). (72) 



D,{z) 

Alternatively, several model calculations of the Collins function have been performed and could be used to obtain 
an estimate of the weighted asymmetry, rather than of an upper bound on it. A review of models has been 
presented in Ref . [40| . The first model was given by Collins [i| and employed in Ref. [4l| to estimate the Collins 
effect asymmetry in SIDIS, which led to the conclusion that H^ '{zi)/Di{zi) = 0{1). Here we will use this 
result and not worry about the sign of the Collins function. Also, we will use that the average value of F{y) is 
approximately 0.5. 

In order to get an estimate of the true asymmetry without artificial enhancement of the weight, one should 
compare Eq. ([S5|) with the weighted cross section {Q\/AM'^), rather than with (1). From Eqs. (jHSj) and (f57| one 
obtains: 

/^\ _ 1 ( Dt''\^,) Dt^\z2) \ 

WpJ-2\ D,{z,) + D,{Z2) )^^^- ^^^^ 

This we estimate by using that 1?^ V-^i = {k'j?){z)/{2z'^NP). Ref. [i2| presented a fit to LEP data to find that 
the average transverse momentum squared of pions inside a jet can be parameterized well as 

{k'^?) 5 (z) = 0.61 z°-2^(l - z)0-20 GeV/c. (74) 

This leads at zi = Z2 = 1/2 (where the average transverse momentum squared is maximal approx- 
imately) to {Q'^/AAI'^^ Ri 20 (1) and hence to an estimate at the few percent level for the ratio 
((Qt/4-M^) cos2(/)i) / ((5|,/4Af^). Of course, this should not be viewed as more than a crude estimate. 

Note that H^ ^ '{zi)/Di{zi) — 0(1) is consistent with the bound in Eq. (|72p . if one uses the maximum of 

(fc^2)2(z) « 0.44 (around z = 1/2) for z(|fcj.|)(z), which leads to i/i^^^^(zi)/i:>i(zi) < 3. 

Another model study [33] has also obtained a prediction for the twice-weighted asymmetry (cf. Eqs. ([S7| and 
(|68l) with Ml = M2) 

,2 oA \m - J(fPh±Pl^cos2(j)id^a 



{P^^COS2(j)i){e2,Zi,Z2) 



2B{y) H^^'^Hf^''' , , 

\y) 1 1 ^75) 



A{y) [B^pf ^d'-^^^D^ 
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in e+e^ ^ hi /12 X. The calculation is based on the Manohar-Georgi model [43] and reproduces the unpolarized 
fragmentation function reasonably well. It leads for fixed Z2 bins to an asymmetry that is almost linearly rising 
as function of zi. For 0.5 < Z2 < 0.8 and zi ~ 0.5 the weighted asymmetry is 3-4%. This is similar in size to 
the estimate obtained above. This indicates that measuring the weighted asymmetry should be feasible for the 
present-day high luminosity electron-positron scattering experiments. 

V. UNIVERSALITY OF THE COLLINS EFFECT 

Applying the results for the Collins function from BELLE in the Collins asymmetry in semi-inclusive DIS, 
assumes the Collins function is universal, i.e. that it is the same for all processes in which it occurs. This assumption 
is not as obvious as it may seem. For T-odd distribution functions it has been shown [4J, |4^ that they are process 
dependent. This follows from their gauge invariant definition as matrix elements of operators involving path- 
ordered exponentials that arc nonlocal off the lightcone. For T-odd fragmentation functions, such as the Collins 
fragmentation function, a similar conclusion was drawn (3l|,|46|. However, for the particular case of e'^e" — > hih2X 
and SIDIS, it has been argued usin g a model, but later also with more general arguments, that the Collins 
fragmentation functions are identical [43, |4^. Recently this conclusion was extended to the process pp —> h jet X 
[43] and argued to be a generic and model- independent result. For the use of the Collins fragmentation functions in 
hadronic collisions it is essential though that factorization holds, which is not established however [SO. ISll |53. l53| . 
A further indication for the universality of the Collins function comes from a spectral analysis of the fragmentation 
correlation function within a spectator model calculation |54| . 

Assuming the universality of the Collins fragmentation function for the processes e+e" -^ hih2X and SIDIS, 
a simultaneous fit to the Collins effect asymmetry data has been performed and a first extraction of transversity 
was obtained [3| . This demonstrates the feasibility of using the Collins effect to access transversity and why it is 
worth doing the Collins effect asymmetry measurement at BELLE. For this reason it becomes also important to 
study other potential contributions to the asymmetry that arises from the Collins effect. This we will do in the 
following sections. 

VI. HIGHER TWIST 

In this section we will study the terms that arise when going beyond leading order in 1/Q. Insertion of the 
leading order parameterization Eq. (|52p of A in the calculation of the diagram shown in Fig. [3] also produces l/Q 
contributions. Such 1/Q contribution can already be generated by simply transforming to a different frame. This 
contribution is not electromagnetically gauge invariant and the full calculation at order 1/Q (6| requires first of 
all, that the correlation function A is parameterized further to include higher twist fragmentation functions [33| 



II^UA,^r-^^E.,n^^^.H.^.ntn'^,. 



M'.kr) = ^lD,^ + H,^ "-^!:gp- + E1 + D^^ + Ho,,ntn\ , (76) 



and second, inclusion of the diagrams in Fig. |4] [2l|. These four diagrams involve one gluon which connects to one 
of the two soft hadronic matrix elements. 

Hence, up to and including order 1/Q the quark fragmentation is described with help of two types of correlation 
functions: the quark correlation function A(Pi;fc) discussed before and the quark-gluon correlation function 
A^(Pi;fc,fci) P: 

AS,^.(Fi;ft,A:O = E/(^(0i^''''^^'"^'"'^W^(^)5^T(y)lA;X)(Pi;X|^^.(^^ (77) 

where A;, ki are the quark momenta and again inclusion of path-ordered exponentials and an averaging over color 
indices are understood. Note that the definition of AJ| includes one power of the strong coupling constant g and 

In a calculation up to subleading order, we only encounter the partly integrated correlation functions 
Jdk^A{Pi;k) and J dk^d'^ki A'^{Pi;k,ki). This allows to express the quark-gluon correlation functions in 
terms of the quark correlation functions with help of the classical equations of motion (e.o.m.) [5^. In the 
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./\/\y 




FIG. 4: Diagrams contributing to e^e annihilation at order 1/Q. 



subleading terms of the cross section one encounters functions indicated with a tilde {E, H, . . . ), which differ 
from the corresponding twist-3 functions {H , E, . . . ) by a twist-2 part, namely 



D-L ^ zDi + D^, 



(78) 
(79) 



H 



"T ^TJ± 



Ml 



2 zHt + H. 



(80) 



For more details we refer to ld|. 

The five diagrams lead to the following expression for the hadron tensor up to and including order 1/Q: 

W"' = 3 / dp-dk+(fpj,d?kj, 6^{pt + kj,- Qt) \ Tr (A{p)"/''A{k)Y) 



^«t/ 



-Tr l^A^ip)rA{kh^-^r ] Tr ^7oA^'(p)7oj r^l^Aik)^ 
+ Tr (A{p)r (loA'XHk)^,) 7''^7.) + Tr (A{p)^^-^rA'l{k)Y 



p+ ,fc"" 



(81) 



The terms with fi± arise from the fermion propagators in the hard part neglecting contributions that will appear 
suppressed by powers of Q^, 



9i + m 



(g+-pi+)7 7" 



(q-pi)2-m2 2{q+-pi + )q- 2q- Q^' 



(82) 






(fci -q )7^ 



q-)q^ 



-2g^ 



QV2' 
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(83) 



As mentioned, one can always integrate out one of the momenta of Aj^(fc, ki) or A^(p,pi) and apply the e.o.m. 
immediately. The quantities A'^{k) and 7oA^ (/c)7o arise from integrating out the second and first argument of 
A^(fc, fci), respectively: 

d^k,A^,Jk,k,) 



I 



d^k,A'X,Jk,,k) 



E 

X 

E 

X 



d^x 






(84) 



e''=-(0|V^.(x)|Pi;X)(Pi;X|5A^(0)^,(0)|0) = (7oA^S),,(fc), (85) 



and similarly for Aj^{p) and 70 A^ {p)lo- 

To obtain the expressions for the symmetric and antisymmetric parts of the hadron tensor we expand all vectors 
in A, A, A^ and A^ in the perpendicular basis (t, z and _L directions). In particular, we reexpress the momenta 
krp and p^ in terms of their perpendicular parts and a part along t and z. For this we need 



Uf up V 

9t ^9±9tp 



T- (fi^'- M it^\h^ 



Q 



it' 



We will refer to these perpendicular projections as e.g. k± rather than /c^^. Thus 

fcj^ = .9^ kfi, — krp 



^^ ■ ''^ (t> + z'^), 



Q 



(86) 



(87) 



and similarly for p± . We note that for these four vectors the two-component perpendicular parts are the same as 
the two-component transverse parts, i.e., k± = krp, etc. The full expression for the hadron tensor is then 



w- 



12ziZ2 / d^krp d^prp 5\pp + kp,-qrp)\ - g^'i^DiD 



— kS^p"} + g^^ fe_L -P^L „±T7^ 



z^f^k"} 



Q 



Zl Ml ' Z2 



-2- 



M1M2 

r.{p U} 

z p_^ 



Q 



-HtH, 



Zl M2 Zl 



The cross section at leading order in a^ but including twist-3 contributions becomes 



dcr(e+e -^ hih2X) 
dzidz2dfld'^qp 



%zlzl\A{y)J^[DiDi] 



+B{y) cos2</)i T 
-C(v)D{y) cos 01 



\ A ft' rZrp ft ' P'T' rC'Ti ' P'T' 

( Ml 



Q 

M2 
' Q 



T 



hkrpD^ — 
Di 



T 



Ml Zl 

hprp 



M1M2 
M2 



Di 

M2 Z2 



Q 

Ml 



T 



hkrp ^H 

Z2 



Ml 



T 



hpp H—± 

M2 Zl 1 



(89) 



This is the result in Eq. (|53p plus an additional cos0 asymmetry of order M/Q. The above expression is given 
in the GJ frame. When transforming to the CS frame the cos 2(f) asymmetry remains unchanged up to terms of 
order MQt/Q^, as can be seen from Eq. ([4T|) . 

The function E only contributes in the case of polarized electrons at 1/Q and for unpolarized electrons at the 
1/Q^ level. The latter would be the analogue of the contribution considered by Jaffe & Ji [5y| for the Drell-Yan 
process after the replacement S'^{prp + kp — qp) —^ S'^{qrp). 
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We have neglected dynamic twist-4 effects, which are order Af^/Q^ corrections. At BELLE these effects are 
expected to be at most at the percent level, but nevertheless could be relevant to include. To inve stig ate this 
further theoretically requires an extensive study due to the many possible sources of l/Q^ effects |28l. |29|. ISa. |57| . 
which thus far has not been undertaken. One reason for this is that factorization has not been considered yet in 
this case (actually not even for the twist-3 case thus far). One particular type of higher twist effect was studied by 
Berger [22,] for e+e^ -^ ttX, which is mostly relevant at large values of z/i, towards the exclusive limit. This is like 
the higher twist contributions to the Drell-Yan azimuthal asymmetries studied in Refs. [5a, [S^, which contribute 
mainly at large x. However, it remains to be seen whether this specific Zh dependence holds true for all types of 
higher twist effects. They are generally expected to lead to /x > i'. Therefore, if /z is found to be much smaller 
than V, it would be a strong indication that one is not dealing with higher twist effects. Experimentally one could 
also test whether such effects are relevant by allowing in the fits to asymmetries for additional 1/Q^ dependence, 
as has been done in studies of DIS data [60, [fill) [£3 ■ 

VII. ELECTROWEAK INTERFERENCE EFFECTS 

In the previous sections we have presented the results of the tree-level calculation of inclusive two-hadron 
production in electron-positron annihilation via one photon up to and including order l/Q, where the scale Q is 
defined by the (timelike) photon momentum q (with Q^ = q'^) and given by Q = \/s- The quantity Q is required 
to be much larger than characteristic hadronic scales, a requirement satisfied by the BELLE experiment for which 
-^s ^ 10.5 GeV. Although this energy is well below the Z-boson mass, we will now consider j-Z interference 
effects, which could lead to percent level contributions (©(s/Af^)). 

Only leading order (l/Q)" effects are discussed, since the combination of power corrections of order 1/Q and 
7-Z interference effects is expected to be negligible (permille level). Here we will only focus on tree level, i.e., 
order (as)°- This may be improved upon at a later stage, following e.g. Ref. [63l |. 

To leading order the expression for the hadron tensor, including quarks and antiquarks, is 



3 I dp-dk+(fpj,d?kj,5^{pj, + fey - Qt) Tr {^{p)V^'l\{k)V'')\ + ,_ + ( ^ ^ ^ | , (90) 



where for a photon V^ = ej'^ and for a Z boson V^ = gvl^ + gAlbl^ ■ We have omitted flavor indices and 
summation. The vector and axial-vector couplings to the Z boson are given by: 

gl, = T|-2QJ"sin2 0w, (91) 

9\ = n^ (92) 

where Q^ denotes the charge and r| the weak isospin of particle j (i.e., r| = +1/2 for j = u and r| = —1/2 for 
j = e-,d,s). 

The lepton tensor is given by (neglecting lepton masses and polarization) 

L%{1, 1') = C"^" [21/^ + 21J'^ - Q^g^,] + D'= 2i e^^pj^r , (93) 

where we have defined 

C^^ = l, C^^ - C^^ - e'si,, C^^ ^ g'v' + g'A\ (94) 

D^7^0, i^^^ = i?^^ - -e'g^, D""^ = -2g'yg'^ , (95) 

where e' denotes the coupling of the photon to the leptons in units of the positron charge; gy, g\ denote the 
vector and axial-vector couplings of the Z boson to the leptons, respectively. 

The cross section is obtained from the hadron tensor after contraction with the lepton tensor, here given in the 
perpendicular basis, 

Lt; = C^^Q^ -{l-2y + 2y^)g^^+Ay{l-y)z^^z'^ 

-Ay{l - y) [l^Jl + \ .<") - 2(1 - 2y),/^(l^) z^^l^^ 
+D^=Q' I (1 - 2y) ef - 2i ^y{l - y) Lpe^}^ z"^ . (96) 
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Here only unsuppressed results arise when both indices /u, v are in the perpendicular directions. 

TABLE II: Contractions of the lepton tensor LJf^ with tensor structures appearing in the hadron tensor. 



i;f.«'''7(4Q') 






C^'^h-y + y") 



-C'^y(l-y)|ax||6±|cos(<?i„ + </.6) 
C'^y(l-y)|ax||b±|sin(,j!.„ + 06) 



Using Table |TT] one finds for the leading order unpolarized cross section, taking into account both photon and 
Z-boson contributions 64], 



"£^:.^^r -l?'''^g {"■'"' "i^-^-i 



+ [K^{y) cos(20i) + K^iy) sin(20i)] T 



Zj ft/ ' P'T~' •" ' ni'Y' Pt' ' '^'T' 



Ml A/2 



The functions Kf{y) (as before, a is the flavor index) are defined as 

Kl{y) - A{y) [el + 2eig'yeag1.Xi + c{clx2] - ^ \leig\eag\xi+ c\cl 
Kl{y) - B{y) [e^ + 2eig'ye,g<^Xi + ^0^X2] , 
Kl{y) = B{y) [2eig'yeag\xA , 
which contain the combinations of the couplings 

4 = {9i.'-9\'), 



X2\ 



The propagator factors are given by 



Xi 



X2 



X3 



3 = 29Wa- 

Q2(g2_M|) 



1 



sin2(20H') (Q'-M|)'+r|M|' 

1 c?2 



-TzMz 
Q2 „ Ml 



Xi- 



(97) 

(98) 

(99) 

(100) 

(101) 
(102) 
(103) 

(104) 
(105) 
(106) 



Using Mz = 91.1876 ± 0.0021 GeV « 91.19 GeV and Vz = 2.4952 ± 0.0023 GeV « 2.50 GeV, and using 
s\T? 6w{Mz)\jjg = 0.23113(15) w 0.231 (and assuming very slow running of this quantity [6^), we obtain for 
g = 10.50 GeV 



Xi w 0.019, X2 ~ 0.00036, xs « -0.00052. 



(107) 
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This shows that we can neglect the sin20i asymmetry in Eq. (|97p and also the ZZ contributions in K\ and K^. 
This leads to the following approximations for K"^ in the lepton center of mass frame: 

e2 
^"(y) ~ -f [1.0005(1 + cos2 6i2) + 0.03 COS 6*2] , (108) 

g2 

Kiiy) w -J [1.002(1 + cos2 6(2)+ 0.06 cos 6I2] , (109) 

K^y) « f [1.0004 sin^ 02], (110) 

Kiiy) « ^ [1.001 sin2f?2], (111) 

g2 

Kiiy) ~ -f [0.00003 sin^ 6I2] «0, (112) 

g2 

Kfiy) ~ -f [0.00006 sin^ 6I2] ~ 0. (113) 



In general, the j-Z interference corrections are small (permille level), except for the forward-backward asymmetry 
(^ cos 6*2), which is on the few percent level. The latter agrees with the results of Ref. [60|, taking into account 
that their angle 6 is defined as the angle between the incoming electron and the outgoing quark, leading to an 
additional minus sign for the asymmetry. 

In order to estimate the size of the hadron level forward-backward asymmetry, we will again consider the 
case of 7r+ and tt^ production from only u and d quarks. Furthermore, we assume D"^'^ (z) = Df^'^ (z), 
Df^^ (z) = Di^^ (z) and neglect unfavored fragmentation functions like D'^^'^ (z), etc. If we define 



da ^ [1 + cos^ 02 + ApB cos 62 + ■■], (114) 

then these assumptions lead to ApB ~ 4%. 

In conclusion, ^-Z interference does not lead to any significant additional (/) dependence at -^s ~ 10.5 GcV. It 
only modifies the 9 distribution with a forward-backward asymmetry term of a few percent. 

VIII. JET FRAME ASYMMETRY 

In an unpublished study [63| a transverse spin correlation similar to the cos 2</) in back-to-back jets was ex- 
perimentally investigated to some extent using LEP's DELPHI data. The following angular dependence of the 
differential cross section for correlated hadron production in opposite jets was studied: 

ex H-cos^6'-|-CTTS'sin^6'cos(0 + (/;>')- (115) 



dcos 

Here ctt = (l^gp— kgP)/(kgP + k9p), with Vq and Og the vector and axial- vector couplings of quarks to the Z 
boson, respectively. So in our notation ctt — C2/c°. S is the analyzing power of the asymmetry to be determined. 
The azimuthal angles (p and <j)' are of the leading particles in the two jets with respect to the qq axis in the lepton 
pair center of mass frame. 

Such a cos{4> + 4>') asymmetry differs from the cos20i asymmetry discussed thus far in that now three momenta 
in the final state need to be determined, namely besides two hadron momenta also the jet axis. Hence there are 
two azimuthal angles, (f) and (/)', that need to be measured. Furthermore, as we will show, the analyzing power S 
is a different expression of the Collins functions. We find that it involves moments of the functions H^ and H^ , 
different from the ones in the correlation Eq. ([55]) that we considered thus far. 

Unlike Ref. [67[ we will not only include the Z boson, but also the photon contributions, including photon-Z 
interference terms (except for the interference term analogous to K^{y) sin20i, which was found to be very small). 

To derive a tree level expression for the analyzing power S in terms of fragmentation functions, we start with 
the hadron tensor expressed in the transverse basis 



W,7 = 12ziZ2 f d^kj. d^PT S^pr + fcy - Qr) 



9t ^ij ~ *^T -^Jj 



DiDi 
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_k_^^^+g£kr:P^ A (116) 

M1M2 «j i i j ' y ' 

where in analogy to the lepton tensor Eq. ([M|) we have defined 

(117) 
(118) 
(119) 

In previous sections we have transformed such hadron tensor expressions into the perpendicular basis of the 
GJ frame, after which the contraction with the lepton tensor yields the cross section. Now we will follow Ref. 
[671 . If one would determine the jet-axis, which is identified with the qq axis, or as an approximation to it the 
thrust-axis, then a measurement of the transverse momenta of the leading particles in the two jets compared to 
the jet momentum is a determination of the transverse momenta of the quarks compared to the leading hadrons 
they fragment into. One can then keep the cross section differential in the azimuthal angles of the transverse 
momentum of the quarks, after which the Qrp integration can be safely done (as opposed to the case of the cos 2(/>i 
asymmetry) and it will not average to zero unless one integrates over the azimuthal angles. In this way one will 
arrive at an expression involving the moments 



^77 


-el 


C^z - C^z, - 


e'-g'^, 


no. _ a 
'-'ZZ - Cl , 


^;7 = 


= 0, 


D^yz = D% = 


-^"9% 


^ZZ — "-3 ' 


^7% = 


2 


^-yZ - ^Z-y - 


e'-g'^, 


^ZZ — ^2 • 



fN(^,)= d\k 



|2 



M, 



F{z,,\kT\^), (120) 



for n = and n = 1. The latter is often referred to as the "half-moment" and also written as F^-^/^\ 

To make the transformation from the frame in which Pi and P2 are collinear (for which we employed the 
transverse basis) into the lepton pair center of mass frame where the qq axis defines the z axis, which we will 
refer to as the jet frame, a new perpendicular basis will be defined (the new perpendicular directions will also be 
indicated by _L). We will choose: 

i^ - f , (121) 



ki" -pi" 



(122) 



We find that 



9T - 5r - ^ [p't'K' + 4""-') - ^pTK'- (123) 



Hence, Pi± = —z\ kx and P2± — ~Z2Pt, up to 1/Q^ corrections. 

Since the z axis is defined differently now, the contraction of the hadron tensor with the lepton tensor yields 
somewhat different expressions. We give the relevant contractions in Table ITTTl The functions A, B and C are the 
same functions oi y — P2 ■ 1/ P2 ■ q as before. If one expresses the lepton momentum / in terms of t and z, i.e., 

; = Qi/2+^Q'^l'i-l\z + l^, one finds that A(y) = l/2-B(y) and C(y) == -^l - A B{y), With B{y) =l\/Q^. 

Using Table [mi we obtain in leading order in 1/Q and as the following expression for the cross section differential 
in (p and 0', in case of unpolarized final state hadrons: 

2 f 

^2) 

>2)|. (124) 

This is to be compared with Eq. (|53p . For simplicity of comparison we will now consider the one-flavor case 
(w-quark dominance). Hence, we find for the analyzing power S of Eq. (jllSp 

Ht^%i)H-^%2) _^ 



da{e+e- -> /ii/i2^) 
dzidz2dfld(f> d(j)' 






+K^{y) cos(<^ -1- cj)') Ht^^^^{zi)H^^^^^ 
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TABLE III: Contractions of the lepton tensor L^j^ with tensor structures appearing in the hadron tensor, for the jet frame. 



iM'^w''7(4Q') 



-3r 4 A{y) 

k^^p^r,} + (feT ■ Pt) 9t "C'l B{y) |Pix | |P2± | COs((/> + (?i')/(^l22) 

i^r 4C(2/)/2 

This can be compared to the expression for the weighted cos20i asymmetry (cf. Eq. ((M|) ): 



where 

i?l(zi)i?l(z2) 

If one now simply assumes a Gaussian fe^-dependence of the functions, i.e., 

H:t{z,k'T^) = H^{z)R^ exp(-i?2fc^)/(7rz2) (128) 

and similarly for Di{z, k'rp^), and equal masses and widths, then one finds that 

o_ ^ HHzi)Hi{z2) 

4R^M^ Di{zi)D,iz2) ^ ' 



and 



4i?4M4 Di{zi)Di{z2) 
Note that in these expressions only the Gaussian width of H^ appears. Therefore, one finds: 



(130) 



For a pion we assume RM k, 0.5, such that S" « 45/ tt. As 5' was expected to be on the few percent level, we 
also expect the analyzing power of the cos((/) + 0') asymmetry in the jet frame to be of that magnitude. 

If one considers Gaussian transverse momentum dependence one does not need to resort to weighting to deal with 
the convolutions. One can consider simply the Ql^-integrated cos 201 asymmetry (numerator and denominator 
integrated separately). One arrives at: 



S" = 



J dq^J^ (2h-krp hpj, — kj,pj, J H^H^ 



IVP J dq^T [DiDi] 
Ht{zi)Hi{z2) 2 



2i?2M2 Di{zi)Di{z2) TT 



S, (132) 



which is somewhat smaller than the analyzing power S of the cos{(j) + 0') asymmetry. All this indicates that 
changing from the GJ frame to the jet frame does not modify the analyzing power of the Collins effect asymmetry 
much. 
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IX. SCALE DEPENDENCE OF THE COLLINS EFFECT ASYMMETRY 

Thus far we have only discussed the azimuthal dependences that arise at tree level. This means that the 
expressions are valid in the region where the observed transverse momentum Qt is small compared to the hard 
scale Q, applicable only in the region of intrinsic transverse momentum. After the extraction of the Collins function 
at BELLE, which is performed at one particular Q^ value, the extracted functions will be used in asymmetries at 
different energies, for instance at lower Q^ for SIDIS or higher Q^ for a comparison to LEPl data. For this one 
would need to know the scale dependence of the Collins function and asymmetry. However, the study of the scale 
dependence of expressions involving TMDs is not far developed. Here we will follow (and improve) the discussion 
of this issue given in Ref . [Ml . 

It is clear that coUinear factorization is not the right framework to address the question of scale dependence of 
the Collins effect asymmetry. Rather, for small Qt the so-called Collins-Soper (CS) factorization theorem J38J is 
of relevance (see also Ref. [69| for a discussion of factorization and transverse momentum in two-hadron inclusive 
e^e~ annihilation). It applies to the differential cross section: 

dzidz2dVt(Pqj, 

Collins & Soper [B^ proved their factorization theorem for the process of interest here: e^e~ — > hi ft,2 X, but 
without inclusion of quark spin effects. A similar factorization for SIDIS and Drell-Yan, including polarization, 
was discussed more recently by Ji, Ma & Yuan [70], with some small differences w.r.t. Ref. [68|. 

In this section we will study the scale dependence of the Collins effect asymmetry taking CS factorization as 
our starting point. It allows us to obtain the dominant Q^ dependence of the asymmetry at small Q^, but the 
analysis also extends the range of applicability of the asymmetry results from the region of intrinsic transverse 
momentum to the region of moderate Qrp values (still under the restriction that Q\ <C Q^). The consequences 
discussed below affect all transverse momentum dependent azimuthal spin asymmetries. It turns out that such 
azimuthal asymmetries generally suffer from Sudakov suppression with increasing Q^ in the region where the 
transverse momentum Qt is much smaller than the large energy scale Q. This Sudakov suppression stems from 
soft gluon radiation. In Ref. ]V{ Collins remarks that Sudakov factors will have the effect of diluting the (Collins 
effect) single spin asymmetry in semi-inclusive DIS, due to broadening of the transverse momentum distribution 
by soft gluon emission. Here we will study this effect in a quantitative way for the cos 2(/) Collins-effect asymmetry 
in e+e~ -^ /11/12X. It shows that tree level estimates of such asymmetries tend to yield overestimates and 
increasingly so with rising energy. 

A. Collins-Soper factorization 

In the CS formalism the differential cross section at small Q^/Q^ is written as 



where 

Wib,Q;zi,Z2) = Y,bl{ziMllh,as{l/h))Y,D\{z2,b-\/h,as{l/b)) 

a h 

xe-SC-'Q) Hab (Q; as{Q)) U{b; 1/b, a.(l/6)). (135) 

Here Di {z, b) is the Fourier transform of the transverse momentum dependent fragmentation function Di {z, z'^krp); 
^-s(b,Q) jg ^j^g Sudakov form factor; H is the partonic hard scattering part; and, U is called the soft factor. The 
fragmentation functions and the soft factor are taken at the scale fi — 1/b. Furthermore, note that there are no 
integrals over momentum fractions, those appear in the large Qt or equivalently, small-& limit only. Hence, the 
Zi in Eq. (|135p are the observed momentum fractions. 

The Sudakov form factor arises due to summation of soft gluon contributions. This is in contrast to more 
inclusive cross sections for which there is often a cancellation of such soft gluon contributions. At values b^ = 
b <C 1/A^, the Sudakov form factor is perturbatively calculable and of the form 



S{b,Q)^ f %- 



Q 



2 






(136) 
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where bo = 2exp(— 7^;) w 1.123 (here we use the usual CS factorization constants Ci = bo, C2 = 1)- One can 
expand the functions A and B in Us/n and the first few coefficients are known, cf. e.g. [7l|, 122]. At the leading 
logarithm (LL) level one needs to take into account only the first term in the expansion of A: A^^^^ = Cp/'n', which 
is the same for unpolarized and polarized scattering. It leads to the double leading logarithmic approximation 
(DLLA) result ^i 



S{b,Q)^CF^log'{b^Q') 



Including the running of as leads to the expression [7^ 



Sib,Q) 



16 



33 - 2nf 



log 



2n2 



b'Q 



bl 1 '^°nA2 



Q' 



log 



iog(b^gVb§) 

log(QVA2) 



(137) 



(138) 



We will take for the number of flavors Uf — 5, since we are interested in energies just below the T(4S'), which is 
above the bb threshold. Furthermore, we take Kqcd — 200 MeV. 

Using only the perturbative expression for the Sudakov factor in the cross section expression (jl34p is valid for Q^ 
very large, when the restriction 6^ <C 1/A^ is justified. If also b^ > \/ A? contributions are important (/x^ < A^), for 
example at small Qt, then one needs to include a nonperturbative Sudakov factor. This can be done for instance 
via the introduction of a 6-regulator [75| : 6^6*= 6/ ■\/l + 6^/6^^, such that 6, stays always smaller than 6max- 
Usually 6max = 0.5 GcV^^, such that as(6o/fe* J^ 2) <0.3. For the function T4^(6*) a perturbative expression for 
the Sudakov factor can thus be used. 

One can rewrite W{b) as: 



W(b) = W{b^)e-^"''^'^\ 
In general the nonperturbative Sudakov factor is of the form [75J 

SNp{b, Q/Qo) = giizi,b) + 52(^2, fo) + ln(QVQo)53(6), 



(139) 



(140) 



where Qo ~ 1/^max is the lowest scale for which one trusts perturbation theory. The gt functions {i = 1, 2, 3) are 
not calculable in perturbation theory and need to be fitted to data. In fact, for the Drell-Yan process they have 
been shown to be necessary to include in order to describe available data [73|- Note that Snp is in general Q^ 
dependent, unlike what one may at first thought expect for a nonperturbative quantity. 

Due to the lack of knowledge of the nonperturbative Sudakov factor for the case of interest here, the quantitative 
results obtained below about the size and Q^ dependence of the suppression should not be taken too literally. 
Nevertheless, one can draw conclusions about what determines the Q^ dependence of the transverse momentum 
distribution of the asymmetries and about the size of the suppression for generic nonperturbative Sudakov factors. 

Concerning the other factors in the CS factorization expression, it should be noted that the one-loop expression 
for H and U were obtained in Ref. [68]. At the scale fi == l/b the 6-dependence of U disappears at order as in 
the MS scheme. 

Due to the integration over 6, the CS factorization expression explicitly requires knowledge of the scale de- 
pendence of the 6-dependent fragmentation functions. In Ref. [6J] this dependence was assumed to be weak 
(logarithmic) and neglected, but here we will include it and use that it is determined by the quark field renormal- 
ization. As already mentioned by Collins & Soper [6g] in some cases it may be more convenient to take instead 
of the varying scale fj, = l/b, a fixed scale Mo ^ 1 GeV. That scale was called fiL by Ji et al. [70|, UM ^'^'^ ^^''^ 
we will call it Qo- Using the renormalization group equations given in Ref. |68i] . the CS factorization expression 
P35p can be rewritten as: 



W{h, Q; zi,Z2) = J2 Dl{zi,b- Qo, «.(Qo)) J2 ^1(^2, b; Qo, a,(Qo)) 

a b 



(141) 



An alternative method has been put forward in Ref. (Tg 
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This has the advantage that the ^-dependent fragmentation functions are always considered at the same scale Qq 
when integrating over b. 

We first consider the above expression in the perturbative regime: b < l/Qo- The one-loop expressions for H, 
U and S are 

HabiQ;asiQ)) (X Si,ael{l + a,{Q^)Fi+0{al)), (142) 

U{b;Qo,as{Qo)) = 1 - ^^^^Cp {log Qlb^ + F2) + 0{al), (143) 



S{b,Q,Qo) = Cf I 
Jc 



d^^ asifi) 



Ql M' 



Q 



2 



2l2 



\og^+\ogQ'ob 



(144) 



Here Fi denote renormalization-scheme-dependent finite terms which will be neglected in the following. Following 
Ref. [73| , one might also consider dropping in the expression between brackets in Eq. (1144p the second term with 
respect to the first one in order to arrive at the double leading log expression. This is because log Q^b^ is in general 
not a large log, since the small-&^ (~ 1/Q^) region contributes little for small Qt- In that approximation S is 
not a function of b anymore. It leads to the double log term in S" ex log^ Q^ IQ\- This means there is no coupling 
between Q^ and Q\ (ignoring for simplicity the Q^ dependence of SNp{b, Q/Qo)), such that the Sudakov factor 
exp(— S') will cancel in asymmetry ratios. As a consequence, azimuthal asymmetries will be Q^-independent in 
this approximation. 

However, the logQo^^ term in S is important to keep in order to obtain the leading Q^-dependence of azimuthal 
asymmetries. If one keeps this term, then one obtains (ignoring the running of as for the moment): 

S{b. Q, Qo) = Cf^ [V Q'b' - V Qlb'] . (145) 

In this expression the second term can safely be ignored, if one is interested in the leading Q^-dependence of 
asymmetries only. The average b^ probed is of order 1/Qt, such that in the small-Q^. region (of order Qq): 
log Q^b'^ ^ log Qob^. The result then coincides with the well-known DLLA expression ()137p . Upon dropping 
all other small logarithms and finite terms, this leads to: 

W{b,Q;z,,Z2) (X Y.elD-,{z,,b;Qo,asiQo))Df{z2,b;Qo,asiQo)) 

a 

X exp f-Cp^^ log" Q^bA . (146) 

This expression, which we will refer to below as the DLLA of W, should capture the dominant Q^-dependence of 
the differential cross section, for Q^ 3> Qq and bearing in mind that this becomes a worse approximation as Q^ 
increases. In the numerical results presented below we will include also the subleading logarithms explicitly given 
above and the one- loop running of a^. 

B. Numerical study of the Q^ dependence of the Collins effect asymmetry 

In order to study the expression Eq. (|134p with W of Eq. (|14ip . we will assume again a Gaussian transverse 
momentum dependence for Di{z, z'^krp): 

Di{z, z'^kl) = Di{z) Rl eM-RlkD/TTz'' = Di{z) g{\kj,\-Ru)/z'', (147) 

such that the Fourier transform is 

b^{z, b") = D,{z) exp (-^) lz\ (148) 

We do not need to worry about any scale dependence of the width _R„, since by construction the fragmentation 
functions are always at scale Qo- Here we will not aim to connect to the large-Qy behavior of the cross section, 
where the Gaussian behavior is certainly not correct anymore. This fact and the assumed factorization of b and z 
dependence implies however that it is not guaranteed that the function Di{z) that appears here is indeed exactly 
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equal to the integrated fragmentation function Di{z; Qq) at the scale Qg. Therefore, the safest thing may be to 
construct ratios of asymmetries where the functions -Di(-z) drop out, cf. Sec. IXIII After these words of caution we 
will make some further pragmatic assumptions in order to be able to numerically investigate the Q^ dependence 
of the Collins effect asymmetry. 

The 6-dependent piece of the fragmentation function Di(z, 5^), taken to be a Gaussian here, can be viewed as 
the Q^-independent part of Snp, because they are simply indistinguishable (ignoring the subtlety of using b here 
instead of 6*). Any remaining z dependence that does not factorize can in principle also be included in Snp- 
We further assume that Smp is flavor and spin independent, i.e. it is the same for all fragmentation functions. 
The assumption of spin independence will get better as Q^ becomes larger. Refinements can be included at a 
later stage if the accuracy of the data demands it. Here we will take for the nonperturbative Sudakov factor 
the parameterization by Ladinsky fc Yuan, which was fitted to the transverse momentum distribution of W/Z 
production in pp (pp) scattering [79|. We will use it with the additional simplifying choice of a;ia;2 = 10~^: 

SNp{b, Q/Qo) = gib' + g2h' In (^\ , (149) 

with gi = 0.11GcV2,g2 = 0.58GeV2,Oo = 1.6 GeV and &,„ax = O.SGeV^i. Note that for Q = 10 GeV the 
Q^-independent part is negligible, which justifies ignoring the above-mentioned subtlety. 



Although quantitatively the results will depend considerably on this choice of Snp (cf. Ref. [6j| for a detailed 
discussion), the Q^ dependence of the asymmetry turns out not to be very sensitive to it. Nevertheless, our results 
underline the importance of a good determination of Snp for the proper extraction of the Collins function, when 
going beyond tree level, see Sec. IIXDI 

The nonperturbative Sudakov factor allows us to integrate W safely over the large b region, but one also has 
to deal with the region of very small b. S{b, Q, Qq) at very small b^ < 1/Q' requires regularization in order to 
ensure the correct limiting behavior. At small Qt this issue becomes less important as Q^ increases. Usually S 
is regulated by the replacement [73| : 

log'(Q262) ^ log'iQH'' + 1). (150) 

In summary, we consider the following pragmatic simplifying steps: we assume a factorized Gaussian dependence 
of the fragmentation functions; we ignore flavor and spin dependence of this Gaussian dependence; we take a 
generic nonperturbative Sudakov factor known from pp scattering; we do not worry about matching to high Qt] 
and, we drop the finite terms F,;. All this leads to a manageable expression for W , which in DLLA becomes: 



W^(6,Q;zi,Z2) ex Y^elDUz^) DHz^) 



X exp (-Cp^^ log2 {QHI + l)\ exp (-5^p(6, Q/Qq)) , (151) 

But as said, below we will also include the single logarithms and the running of a^, and doing so leads to 
qualitatively similar, but quantitatively different results compared to earlier results of [64|, where ^ = 1/6 was 
used. 

The final step to be taken is to include the Collins fragmentation function. This can be done via the replacement 
of D(z, b^) —>■ A(z, 6), where for unpolarized hadron production: 

M^,b)^^S^D,(z,b-)^+[^^H-iz,b-)^ I?}, (152) 

which is the Fourier transform of Eq. ([5^ (we have included a factor P^ into the definition of A). Since the 
second term is b-odd, it leads to a different Q^-dependence than the first term. This leads to the Q^ dependence 
of the asymmetry. 

A model for the transverse momentum dependence of the function H^ is needed. As for Di we will simply 
assume a Gaussian form: H^°'{z, z'^k^) = H^°'{z) Q{\kj,\; R)/ z"^. Strictly speaking, the radius R should be taken 
larger than i?„ of the unpolarized function Di, such as to satisfy the bound [33] 

Ifc^l \Ht(z, IfcrDl < Mh Di{z, Ifc^l), (153) 

for all \krp\. But since we will include this (in principle spin dependent) Gaussian dependence into Snp, we will 
approximate R by i?„, because the (spin independent) Q^ dependent part of Snp dominates anyway. Note that 
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equating R by i?„ is not allowed in tree level treatments of the asymmetry, as no fall-off with Qt will be obtained 
otherwise. 

As a further simplification, we will assume that the fragmentation functions for both hadrons are Gaussians of 
equal width, i.e. we take Rui = Ru2 — Ru and i?i = R2 = R = Ru- Also we take Mi = M2 = M. All these 
simplifications can be easily undone when needed. 

With all these steps in place, we can write the Collins effect asymmetry in the CS formalism in a form similar to 
that obtained at tree level in Sec. IIVBI If we define the asymmetry A{qrp) as the analyzing power of the cos 2(j>i 
asymmetry : 

d.MJ,l ^ °-{l + cos2^iA(g,)}, (154) 

then at tree level we obtained 

^^^^ QlM,M2Y.bKl{v)T[D,D,] ■ ^ ' 

Here we taken the fiavor summation out of the definition of JF[. . .] and replaced e^Aiy) -^ Kf(y) and e^^B{y) -^ 
K^{y) compared to Sec. IIVBI in order to include 'y-Z interference effects, although we do ignore the sin2(^i 
Collins effect asymmetry that was estimated to be small in Sec. lIVDl We have also multiplied the asymmetry by 
a trivial factor Q^/Q^ in order to be able to replace h — > Qrp. This might seem problematic at Qt = 0; but that 
turns out not to be a problem as the asymmetry has a kinematic zero at Qt — 0, because h cannot be defined in 
that case. 
At tree level 

T[DD]= I (fpj.(fkj.5'^{p,^ + kj,-qj,)D°'{zx,zlp%)D''{z2,zlkl). (156) 

In order to apply the CS factorization expression we can replace in Eq. (|156p [64| 

5^{p^ + k^-q^)^J^ ,-MPr+K-qr) e-S U, (157) 

leading to (suppressing the flavor indices and the arguments of 5* and U) 

e'b.q^e'^UD,{z,,b^)^i{z2,b') 

= — / dbb JoibQr) e-'^ U D{zi,b) D{z2,b). (158) 



The numerator in Eq. (jl55p . 



T [(2qrp-pj, qj,-krp — q^prp-krp^ H^H^] = j 



<Pb 



ib-q^-s, 



(2^)2 
X I (fprp (fkj. {2qj.-prp qrpkj. - qj^Pr-k^) c-'^-'^Pt+^t) H^{zj^^z^p'^)H^{z2, zjk^), (159) 

cannot be treated exactly like the denominator. But for Gaussian transverse momentum dependence of the 
fragmentation functions the transverse momentum integrals can be performed. One finds 

/ (fpj. (fkj. {2qT-PT qr-kr - QxPT-kx) e^'^'^PT+kr) g{p'^;R) g{k'^;R) = 

- ^ [2(Qt • b)' q'rb'] exp (^-^) , (160) 

which after application to Eq. p59p yields 

J'[{2qr-PTqT-kT-qTPT-kT) HtHi] = j A^e^b-q^ (-^) [2{qT ■ bf - q^b"] 

xe-^U Ht[zi,b^)Tl\{z2,b''). (161) 
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It is important to note that the factor l/i?* stems from the Gaussian width of the functions H^{z,krp) to be 
used in this expression. Therefore, it is not dependent on the scale Q as discussed above aheady. 
Putting everything together Eq. (|155p can be transformed into 



AiQr) 



E, i^f (2/)i?f(zi) Hi-{z, 



AM^R^ Zb ^i(y)^i(2i) Di{z2 



=rr^AiQ^ 



EbKl{y)D\{zi)D\{z2) 



AiQr), 



where 
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We wiU employ this expression in combination with Eqs. p43p . (jl44p and p49p . dropping the finite terms Fi, but 
including the one-loop running of as ■ As the replacement of Eq. ()150p turns out to have only a minor effect in 
this expression, in contrast to the DLLA expression, it will not be included. 

In Fig. [5] the asymmetry factor A{Qt) is given at the scales Q = 10 GeV and Q = 90 GeV, in order to compare 
the results for the BELLE and LEPl scales. The solid curves are obtained with the method explained here, 
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FIG. 5: The asymmetry factor A{Qt) (in units of M^) at Q = 10 GeV and Q — 90 GeV. The solid curves are obtained 
with the method explained in the text; the dashed-dotted curves are from the earlier analysis of Ref. [64| . 



whereas the dashed-dotted curves are the results obtained from the analysis of Ref. [6^, where fi = 1/6 was 
employed and the scale dependence of the fragmentation functions was ignored. One observes a reduction w.r.t. 
earlier results, but the large Sudakov suppression with increasing Q remains equally strong, i.e. approximately 
l/Q. The DLLA result (not displayed) decreases more slowly, but as said it becomes a worse approximation as 
Q^ increases. 

The cos 20 asymmetry has been studied using DELPHI data {y/s = Mz) and the magnitude was found to be 
small Q- This may have several reasons, one of which could be the Sudakov suppression discussed here. In that 
case a small result at LEPl energy does not imply that also at BELLE the asymmetry has to be small. In Ref. 
[SOJ I a comparison is made of the Collins functions extracted from the DELPHI and BELLE data (here it should 
be emphasized that the DELPHI data analysis remains preliminary and does not consider possible systematic 
effects). At higher Q^ the extraction of the Collins function from a tree level expression becomes less accurate, 
as we will discuss in the next subsection. Therefore, it would be interesting to also compare the asymmetries 
directly (rather than the extracted Collins functions). This could serve as a check of the Q^ dependence of the 
asymmetry as a whole and thus of the CS formalism. 

Of course, one also wants to compare to Collins effect asymmetries in SIDIS at lower energies (e.g. for HERMES 
(Q^) = 2.41 GeV^). The question is what to do at smaller values of Q^l The logarithms that are resummed in 
the CS formalism are not that large to begin with and neither is the relevant 6-range (set by Qo and Q). In this 
case one can set Qq = Q, which means S{b,Q,Q) = 1 and SNp{b) Q^-independent, and consider Qt ^ M. A 
reduction to the tree level form occurs up to small logarithmic corrections of order UsiQ"^) \ogQ^ / Q^ . Tree level 
analyses should yield reasonable results in this case. 
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C. Comparison to tree level 



We will compare the above result for Q — lOGeV with the tree level result (of. Eq. ([5^ ). In the tree level 
expression for the asymmetry it is important to keep Gaussians in numerator and denominator different, in order 
to ensure the bound given in Eq. (|153p is satisfied and an asymmetry is obtained that falls off at larger Qt- 
The tree level expressions for A{Qt) and A{Qt) will be denoted by A''^\Qt) and A'^^^Qt)- They are given by 
(ignoring electroweak interference effects for simplicity) 



A'^'>\Qt)^ 
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and 



^(0) (Qt) = exp [-(i?2 - i?2)Q2^/2] M^QIr'^/RI. (165) 

In Fig. [6] we have displayed the comparison oi A{Qt) at Q = lOGeV and the tree level quantity A^^'{Qt) using 
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FIG. 6: The asymmetry factor A(CIt) at Q = 10 GeV (solid curve) and the tree level quantity Ac'^' (Qt) using R^ — 1 GeV~^ 
and E? /R^ — 3/2. Both factors are given in units of M^. 

the values R'^ — 1 GeV^^ and R^ / R^ — 3/2, which were chosen such as to minimize the magnitude of A^'^\Qt)^ 
cf. [6J| for further discussion. We conclude that inclusion of Sudakov factors has the effect of suppressing the 
tree level result roughly by a factor of 5, whereas for Q = 90 GeV it is more than an order of magnitude. Tree 
level extractions of the Collins function at large Q^ therefore can significantly underestimate its actual magnitude 
(roughly by the square- root of the Sudakov suppression factor of the asymmetry) . It is important to keep this in 
mind when comparing predictions or fits of transverse momentum dependent azimuthal spin asymmetries based 
on tree level expressions applied at different energies. 

The above also shows that upon including Sudakov factors one retrieves parton model or tree level characteristics 
(also noted in Ref. [75|), but with transverse momentum spreads that are significantly larger than would be 
expected from intrinsic transverse momentum (this is supported by the presently available parameterizations of 
Snp in various processes, which usually have Gaussian ^-dependence with widths that increase with Q^, cf. e.g. 
[81]). 



D. Nonperturbative Sudakov factor from BELLE 



Since the previous results depend on the input for the nonperturbative Sudakov factors Snp, which (as a 
function of 21,22) is not determined for the process e+e~ -^ hi h2 X, the numerical conclusions about the size 
and Q^ dependence of the suppression should be viewed as generic, not as precise predictions. Therefore, we would 
like to stress the need for an extraction of the nonperturbative Sudakov factor from the process e~^e~ — *■ hi /12 X. 
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Considering the wealth of data from BELLE this should pose no problem. Here we will give a brief outline of 
how this could be done. 

Thus far the nonperturbative Sudakov factor has only been obtained from older e^e^ —t hi h2 X data 
for the energy-energy correlation function [84I, \Yl,h h J dziZidz2Z2Q^da/dzidz2dQ'^, at various values of Q^ 



•/hi , 

A method based on CSS factorization [8j| was used as discussed in detail in Ref. [75|. CSS factorization 
can be obtained from CS factorization for the^-integrated cross section which receives only contributions from 
unpolarized quarks (cf. the discussion in Ref. [3^ on the relation between CSS and CS factorization for SIDIS). In 
CSS factorization only collinear fragmentation functions appear, which do not affect the energy-energy correlation 
function on account of the momentum sum rule ^^^ / dzzD^ [z) = 1. 

It needs to be mentioned however that the Qt distribution for the Drell-Yan process requires an Snp that is 
dependent on the momentum fractions [77|. It would therefore be better to extract Snp from e^e^ -^ hi h2 X 
data at given zi and 22 , rather than from the energy-energy correlation function. 

To obtain the nonperturbative Sudakov factor the cross section differential in zi, Z2 and Qt needs to be fitted 
(as said, one can integrate over 0i). Either one uses the CSS formalism as explained in Ref. [7a| or one can use 
the CS formalism with a suitable b dependence of the fragmentation functions. Here we will do the latter and 
stay within the pragmatic approach adopted before, employing a Gaussian form (|148p . In this way we arrive at 
(leaving the Qq dependence implicit) 
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This expression together with Eq. ()138p for S{b) and the general parameterization of Snp in Eq. (I140|) can be 
used to obtain a fit of the parameters gi to the data at small Qt (^ Q^, in order to avoid inclusion of 0{Q^/Q^) 
terms that were dropped in Eq. (|134p V Such an extraction would be of more general interest as well, since there is 
a relation between the dominant Q^-dependent part of nonperturbative Sudakov factors in the Drell-Yan process, 
SIDIS and e+e" -^ hih2 X [8II], which would be interesting to test. 

X. RADIATIVE CORRECTIONS 

Thus far we have restricted the discussion to Q^ -^ Q^. In this section we will look at the high-Q^ region, 
where one expects collinear factorization and fixed order perturbation theory to yield a good description of the 
cross section. 

It is well known that gluon bremsstrahlung also leads to angular asymmetries, as the radiating quark or 
antiquark is recoiling away from the original quark-antiquark axis. In this section we will look at this effect at 
order as, i.e. the radiation of an additional gluon into the final state: e^ e^ -^ q q g (we will not include 7-Z 
interference terms; this can be done at a later stage, when required), see Fig. [T] As it turns out the expressions 
are the simplest in the Collins-Soper frame. This frame is the lepton-pair center of mass frame with the z axis 
defined as pointing in the direction that bisects the three-vectors P2 and —Pi and was discussed in Sec. IIII Al 





FIG. 7: Lowest order gluon bremsstrahlung diagrams. 

In the collinear factorization approach contributions to Wt,Wl,Wa and Waa will be generated at order as- 
In the CS frame certain ratios of these structure functions will be independent of the subsequent fragmentation of 
the quark and antiquark. The expressions we will give here are completely analogous to the Drell-Yan expressions 
obtained in Ref. ^]. They are valid for large Qt (^ M), but still restricting to two jet events. Here we will 
not consider the case where one of the hadrons comes from the gluon fragmenting, which becomes relevant when 
Qt '^ Q- For pronounced three jet events it would certainly have to be taken into account. 

At large Qt at least one of the hadron momenta is deviating considerably from the z direction. The hadron 
momenta are again Pi and P2. Upon neglecting Adf /Q"^ corrections, the photon momentum is parameterized as 
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q = Pi/Ci + P2/C2 + qr (cf. Eqs. ©-dH)), such that Q^ = s/(CiC2) - Qt, where 5 = (Pi + Ps)^- It foUows that 
Zi = Q Q^ IQ^ ■ The calculation is done in coUinear factorization, which means that parton transverse momenta 
w.r.t. the produced hadrons are neglected. The parton momenta are defined such that p; = Pi/£,i- 



Following the calculation of [25], but now applied to process e+ e 



+ ^- 



calculate the contribution of the two diagrams in Fig. [7]), one finds in the CS fn 



q q g -^ hi h2 X (that is, one has to 
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where the function K{(i, (2, Q'r/s) is given by 
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This expression shows that only the partonic Wa depends on the lightcone momentum fractions of partons and 
therefore that only the hadronic Wa will depend on the fragmentation functions Di and Di. This is only the 
case in the CoUins-Soper frame. In the Gottfried- Jackson frame all angular dependences will depend on the 
fragmentation functions. 

Using the other standard notation for the angular dependences Eq. ((5^ one finds that the above yields: 



Acs = 



i^cs 



Q2 - iq: 



Qt 



n 



Q2 + |q: 



We see that also in e+e -^ hi h2 X the so-called Lam- Tung relation holds [231124 

1 - Acs - 2j^cs == , 



(169) 
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or equivalently, Wl = 2Waa- Using Eq. (|¥T|l one sees that the relation also holds in the GJ frame and as shown for 
the Drell-Yan process, it is expected that it continues to hold if one includes the gluon fragmentation contribution, 
although in that case both A and z/ will depend on the fragmentation functions (cf. Ref. [85| for a discussion of 
these aspects for the Drell-Yan process, including the effects of resummation) . We also note that the Lam- Tung 
relation is purely a 0{as) result, i.e. it does not apply beyond leading order. Surprisingly, in the Drell-Yan process 
the Lam- Tung relation is known to be violated by much more than the 0(0^) contribution [86|. The distribution 
function analogue of the Collins effect has been suggested as an explanation for this large deviation from the 
NLO pQCD result [83|- It would be very interesting to check experimentally in e+e^ -^ hi h2 X whether the 
Lam- Tung relation is violated as much as in the Drell-Yan process. 

A natural question is how to combine this coUinear factorization fixed order perturbative result, which is valid 
for A'P ^ Q^, with the CS factorization Collins effect result, which is valid for Q^ ^ Q^. This question was 
addressed recently in Ref. '35| for the analogous situation in semi- inclusive DIS. There it was pointed out that 
the two results may simply be added, which gives an expression for the asymmetry that is correct up to power 
suppressed corrections (of order Q\/Q^ in the low Qt region and of order M'^ IQ\. in the high Qt region). It 
is based on the fact that the denominators of the asymmetry expressions in Eqs. ((59|) and (|170p coincide in the 
intermediate Qt region (M^ <C Q\ <$^ Q^) (a fact that is also used in the derivation of the CSS formalism [S^] 
from CS factorization). In addition, one uses the following observations. The above calculation shows that the 
fixed order perturbative expression for u is of order Qt/Q^ when Q^ ^ Q^, hence power suppressed w.r.t. the 
Collins asymmetry expression given in Eq. ([59]) . Conversely, considering the latter in the large Qt (^ M) limit 
beyond tree level (within CS factorization) leads to a result of order M'^/Q^ [39|, which is power suppressed 
w.r.t. V in Eq. (|170p . This implies that upon neglecting power suppressed terms the sum of the two contributions 
reduces to Eq. ([K5|l at low Qt and to Eq. p70p at high Qt- In the intermediate region the results can be added 
because they have a common denominator. For further details we refer to |35| . 

Despite the fact that the asymmetry for all Qt can be described by the sum of the Collins effect expression 
and the fixed order coUinear factorization expression, the latter complicates the extraction of the Collins function 
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and the weighting with Q^, as explained next. Therefore, one wants to construct quantities in which the hard 
gluon radiation component of the cos 20 asymmetry is (largely) absent. One such option is to consider ratios of 
asymmetries where the hard gluon radiation component of the cos2</) asymmetry drops out. The other would be 
to exploit the Lam- Tung relation. Both options will be discussed below. 

XI. WEIGHTED ASYMMETRY BEYOND TREE LEVEL 

In this section we return to the issue of integrating over Qt and weighting with Q^. In the cross section kept 
differential in the transverse momentum Qj, the Collins effect gives rise to a cos 20i asymmetry (both in the GJ 
and CS frames). In principle one can consider this quantity integrated over the length Qt of qrp, but this yields a 
quantity that does not arise in other processes in the same way. Without further assumptions about the Collins 
fragmentation function, such a quantity would only be useful to test whether it has the expected Q^ behavior. 
This formed the main reason for the suggestion to weight the Q^ integration with an additional factor of Q^. 
At tree level this yields an expression involving the first transverse moment of the Collins function. Beyond tree 
level but considering only the low-Qy expression of the CS formalism this remains true and moreover, a quantity 
results that does not suS'er from Sudakov suppression^, as was pointed out in Ref. [6J|. These nice properties 
are spoiled however when the perturbative high-Qr contribution is included. Unfortunately the weighted cross 
section becomes sensitive mainly to the large Qrp contribution (note that we are discussing the asymmetry in the 
weighted cross section, i.e. we integrate numerator and denominator of the Q^-dependent asymmetry separately). 
Because the latter is calculable perturbatively and can simply be added to the Iow-Qt result, it can in principle 
be subtracted from the experimental result before extracting H^ 

One way to subtract it experimentally is to consider ratios of asymmetries where the hard gluon radiation 
component of the cos 20 asymmetry drops out. This is discussed in the next section. The other possibility would 
be to exploit the Lam- Tung relation. Here the idea is to use the structure functions themselves, as opposed to 
the ratios of structure functions A and v. We observe that the quantity 

I d\T^l{W,-2W^^)^^-^^Y.^l Hi^^%z,)H'[^'^\z,) (172) 

does not receive radiative corrections at order a^. This applies in both the GJ and CS frames and will hold true 
even upon inclusion of the gluon fragmentation contribution. 

Once the first transverse moment of the Collins function has been extracted there is another potentially inter- 
esting test of the Collins effect in the CS formalism, which may turn out to be scale independent. This exploits 
the so-called Schafer-Teryaev sum rule [88] 

^ f dz z hI^^\z) ^ 0. (173) 

h '' 

If we define for each hadron type h the function Cu that is a function of Mh and scale /x only, as 

Ch{Mh,^i)= fdzzH^^^\z), (174) 

then one can conclude that for each hadron type this has the same /i dependence (if any) . This is because there 
is no gluon contribution, such that there is no mixing among gluon and quark functions. For each quark flavor 
one has the same evolution properties. Assuming autonomous evolution of the function H^ ^ i^), this implies 
that ratios of C^'s for different hadrons types, e.g. Ctt/Ck, will be scale independent. This could then be tested 
experimentally without the need to know the scale dependence of the Collins asymmetries. Measuring the ratios 
at BELLE, HERMES, LEPl, etc, should then all give the same answer. The only experimental difficulty would 
be to obtain a sufficient coverage in z and Qt to be able to extract the integrated quantities Ch reliably. In Ref. 
[89| an autonomous evolution equation for zH^ (z) has been obtained in the large Nc limit. However, since its 



^ One probes the Sudakov factor at the point 6 = where it vanishes, although that cannot be seen from the perturbative expression 
for S{b) in Eq. I I136I I. which is only valid for b > bo/Q. 
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derivation employed so-called Lorentz invariance relations the correctness of which has afterwards been disputed 
[90| , the assumption of autonomous evolution may turn out to be unwarranted. This issue should first be settled 
before any conclusions can be drawn from the ratios Chi l^h^ ■ But if scale independence of the ratios Chi l^h^ can 
be established, it would offer an experimental consistency check of the extracted Collins fragmentation functions. 

XII. RATIOS OF ASYMMETRIES 

There are several reasons why it may be useful to consider taking ratios of asymmetries for different sets of 
hadrons. These may be experimental reasons such as to cancel systematic effects, but also theoretical reasons 
namely that effects drop out that are independent of the type of hadron considered. 

If the cos 20 asymmetry itself is not that large, the asymmetry in a ratio of two cos 1(^ asymmetries is to good 
approximation also a cos 20 asymmetry: 

dzidz2dnd^q^ ^ ^^-' ' 

As can be seen, effects that are independent of the type of hadrons hi cancel in the difference A''^''^ — A^^^'^ . 
Since fragmentation functions generally do depend on the type of hadron, only effects that are independent of the 
fragmentation functions will cancel. In Sec. |X] it was discussed that in the CoUins-Soper frame the perturbative 
order-as contribution to v (Eq. ()170p ) is independent of the fragmentation functions. This assumes that the 
hadrons arise from quark fragmentation, which should be a good approximation for Qt not too close to Q. 
Hence, under these conditions (in the CS frame for moderate Qt) the order-a^ radiative correction term cancels 
in the ratio Eq. p75p . If in addition ^cs ~ 0, then the conclusion also extends to other frames that are rotations 
of the CS frame, such as the GJ frame or the jet frame discussed in Sec. IVIIII 

Since at moderate Qt the dominant contribution to v from radiative corrections (Eq. (|170p ) is already small 
itself, any uncanceled remainder from higher order corrections or from the gluon fragmentation contribution is 
expected to be negligible. This can be verified via explicit prediction of the radiative correction to the cross 
section, e.g. via a Monte-Carlo study. Due to the cancellation of the dominant radiative correction, the ratio 
method presents a good way to access the Collins effect without taking into account the radiative corrections 
explicitly. The use of different types of hadron pairs, such as tt+tt^, tt^tt^ and tt^tt^, moreover allows one to 
learn about favored versus unfavored fragmentation functions 80] . 

Predictions for the ratio of jet frame cos(0 + 0') asymmetries (cf. Eq. p24p or Ref. 8]) for unlike-sign pion 
pairs to like-sign pairs have been given in Ref. [9l| employing a spectator model. Asymmetry ratio's of the order 
of the experimental data could be obtained, but with considerable uncertainties. Interestingly, kaon asymmetries 
were predicted to be larger than the pion asymmetries. 

XIII. BEAM POLARIZATION 

Thus far we have assumed that no transverse beam polarization is present. However, it is well-known that 
charged particles circulating in a magnetic field become polarized transversely to the beam direction due to 
emission of spin-flipping synchrotron radiation: the Sokolov-Ternov effect [92]. This effect can be significant for 
electrons and positrons due to their small mass. The transverse polarization of particles circulating in a uniform 
magnetic field is the following function of time: 

where 7 — E/ra is the Lorentz factor, p is the bending radius, and a = 5v3/8. More general situations are 
reviewed in Ref. [:93j. The polarization has a strong dependence on the mass in the exponent: the lighter the 
particle, the faster it becomes polarized. 

It is also well-known [63, [OJ, 123, ISa] that beam polarization can affect the angular distribution of produced 
hadrons in e~^e~ —> hX. As explained in '94'| it leads to a cos 20 asymmetry, where the angle is defined w.r.t. 
the beam and beam spin directions. This asymmetry has in common with the Collins effect asymmetry that both 
are transverse spin asymmetries: the former concerns lepton spins, the latter quark spins. Both asymmetries 
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arise due to interference between contributions of ±1 photon helicity states. In the former the magnetic field 
determines the electron and positron spin directions, which then determine the photon polarization state and 
hence, the subsequent decay of the photon into quark-antiquark pairs, which is visible in the angular distribution 
of final state hadrons. In the Collins asymmetry case, the angular distribution of two final state hadrons w.r.t. 
each other filters or selects the quark and anti-quark spin directions, which via the photon are correlated to 
the electron-positron angular distribution. The Collins effect does not contribute to the angular distribution 
of produced hadrons in e+e^ -^ hX, but transverse beam polarization may affect the angular distribution in 
e+e" —^ hi /i2 X. 

Turning to the situation at BELLE: one should learn whether the beams are polarized due to the Sokolov-Ternov 
effect, because it may lead to angular asymmetries that require correcting for. From the beam energies and bending 
radii of the KEK accelerator one can try to estimate the polarization build-up time from Eq. (|176p . Here one 
should keep in mind that the accelerator consists mostly of straight sections. For the low (high) energy ring the 
bending radius is 16.3 m (104.5 m), but the length of bending is 0.915 m (5.86 m), whereas the circumference of 
the rings is a little over 3000 m. Therefore, the polarization build-up time interval is much reduced compared to 
simply using the quoted bending radii in Eq. (jl76p . Moreover, depolarization effects should be very important, 
especially since the interaction point is in the middle of a straight section. Negligible beam polarization at the 
collision point is thus expected. 

The degree of polarization P can be measured experimentally via the process e+e^ -^ /^^M^ for which the 
cross section is [QJ] 

where aL/crr = 4to^/(3^ « is used. Experimental investigations of this cross section can show the presence or 
absence of any significant transverse beam polarization at BELLE. 

XIV. WEAK DECAYS BACKGROUND 

Another test of systematic effects is offered by the process e+e~ -^ t^t^ -^ Tr'^DrTr~i^r (and other such weak 
decays). In this process a cos 20 asymmetry analogous to the Collins effect asymmetry can arise. In this case 
it is calculable within the Standard Model, using the parity violation parameters (usually denoted by a) for the 

tau decays {a-^i = Tl for r -^ it v t)- The self-analyzing parity- violating weak decay of the t's plays a role 
similar to the Collins effect. It correlates the spatial distribution of the produced pions to the spins of the r's. 
The asymmetry does not depend on the r's being transversely polarized on average (in that case the effect would 
have been of order m^/s). 

The analyzing power of the cos 20 asymmetry is not a very simple expression, because it depends on the decay 
kinematics, for details we refer to Refs. [93, llaH- The asymmetry is proportional to aT^+aT^-CxT = ^Ctt, 
where Ctt is given in terms of the axial and vector coupling of the r lepton to the neutral gauge boson: 

which for the Z boson is close to -1-1 and for the ph oton is —1. This change of sign implies a sign difference between 
the asymmetry measured at LEP |lOOl . llOll . Il02l | and what would be obtained at BELLE. The LEP results are 
in good agreement with the standard model prediction, therefore, assuming the same applies to BELLE, it could 
provide a way to check for systematic effects. The same would apply to heavy-quark weak decays. 

XV. SUMMARY 

We have given an overview of the Collins effect asymmetry in the process e+e" -^ hi ft,2 X, where the two final 
state hadrons belong to opposite jets. We restricted the discussion to two jet events and discussed the particular 
situation of BELLE, but many results also are applicable at other e+e^ collider facilities. The Collins effect -a 
spin effect in the quark fragmentation process- gives rise to an azimuthal asymmetry, which in the Gottfried- 
Jackson or CoUins-Soper frames is a cos 2(f) asymmetry. It was our main objective to study to which extent the 
analyzing power of this angular asymmetry is a measure of the Collins effect. 
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At low values of the observed transverse monientum Qt of the photon we employ Collins-Soper factorization, 
which leads to expressions involving parton transverse momentum dependent fragmentation functions. It is within 
this framework that the Collins effect fragmentation function arises. At tree level and leading twist it leads to 
the only contribution to the cos 2(j) asymmetry. Twist-3 effects (order M /Q terms) were shown in the Gottfried- 
Jackson frame to only contribute to /i, i.e. to give rise to a cos (j) asymmetry. Throughout this paper we have 
neglected twist-4 effects that are of order IvP /Q^. They can contribute to the cos 20 asymmetry, but are expected 
to lead to z^ < ^, i.e. a larger cos than cos 20 asymmetry. This property allows to estimate their importance at 
BELLE. 

Beyond tree level, but still at low Qt, the Collins-Soper factorization expression also dictates the scale Q^ 
dependence of the Collins effect asymmetry. This was studied extensively and it was demonstrated that the 
Collins effect asymmetry suffers from considerable Sudakov suppression as Q^ increases. This was already pointed 
out in Ref. [6J], but the analysis presented here improves on that treatment and leads to quantitatively different 
but qualitatively very similar results. We conclude that inclusion of the dominant double and single logarithmic 
terms at order as has the effect of suppressing the tree level result. In the example considered this suppression was 
roughly by a factor of 5 for Q = 10 GeV as relevant for the BELLE experiment. Therefore, tree level extractions 
of the Collins function can significantly underestimate its actual magnitude. It is important to keep this in mind 
when comparing predictions or fits of transverse momentum dependent azimuthal spin asymmetries based on 
tree level expressions but applied at different energies. To properly describe the cross section beyond tree level 
within the Collins-Soper formalism, a good experimental determination of the nonperturbative Sudakov factor is 
required. An explanation of how this can be done at BELLE was given. 

In the Collins-Soper formalism the Collins fragmentation functions enter in a transverse momentum convolution 
expression. Such convolutions do not arise in other processes in the same way, therefore one has to resort to other 
means of extracting information that can be applied in other processes. Often a Gaussian transverse momentum 
dependence of the Collins function is assumed because in that case all integrations can be performed analytically. 
But this introduces a model dependence in the results. Another way suggested is to consider weighted integration 
over the observed transverse momentum Qt- In this way the Q^n-weighted cross section at tree level probes the 
first transverse moment of the Collins fragmentation function, which also arises in other processes identically. 
This would allow for model independent predictions once the Collins functions are extracted from the BELLE 
data. Another advantage of the Q^-weighted cross section is that it does not suffer from the Sudakov suppression 
mentioned above. However, considering the integral over all Qt requires that one also describes the asymmetry 
well at high Qt- This complicates matters considerably. 

For the cos 20 asymmetry at high Qt collinear factorization can be used to describe it. This is most conveniently 
done in the Collins-Soper frame, where the dependence on the fragmentation functions drops out to a large 
extent. The asymmetry arises from hard gluon emission and is not related to the Collins effect. Therefore, 
it becomes important to describe this dependence, that behaves as Q^/Q^, as best as possible in order to 
distinguish it from the Collins effect contribution that at large Qt behaves as M'^ /Q\, (as opposed to Gaussian 
transverse momentum). Due to this specific Q^ dependence of the two effects, it is possible (upon ignoring twist-4 
contributions) to simply add the two contributions (cf. Ref. [3q|). The Collins effect then dominates at low Qt 
and the hard gluon emission at large Qt- Unfortunately this implies that the Ql^-weighted integration is sensitive 
mostly to the latter effect and not to the Collins effect of interest. So one first has to subtract the perturbative 
collinear factorization result at large Qt- This can be done to a large extent automatically by considering ratios of 
asymmetries, but a new way exploiting the analogue of the Lam- Tung relation was pointed out. This relation was 
shown to be violated strongly in the Drell-Yan process w.r.t. the NLO pQCD result. It would be very interesting 
to see whether this also holds in e+e" collisions. 

To get an idea about the size of the Collins effect asymmetry in e+e^ — > hi /12 X, a crude estimate of the Q^- 
weighted Collins asymmetry at tree level was made. It shows the asymmetry can be on the level of a few percent. 
This is in agreement with a model prediction in the literature. These estimates employ tree level expressions, but 
they should remain reasonable estimates beyond tree level (due to the absence of Sudakov suppression), provided 
the dominant high-Qy gluon emission contribution is subtracted. 

A comparison to the Collins asymmetry in the so-called jet frame (defined with help of the qq axis or approx- 
imately by the thrust axis) was given. In this frame the Collins asymmetry becomes a cos(0 -I- 0') asymmetry. 
It was shown to probe so-called half-moments of the Collins fragmentation functions, as opposed to the first 
transverse moments appearing in the (5|,-weighted Collins asymmetry. It is estimated to be similar in size as the 
asymmetry in the GJ frame. 

Other potential contributions to the Collins effect asymmetry have been investigated, such as from 7-Z inter- 
ference and the effect of beam polarization. Both effects should be negligible at BELLE (well below the percent 
level). Elcctroweak contributions can in principle lead to a sin20i asymmetry, but also that was estimated to be 
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very small. Therefore, we conclude that 7-Z interference does not lead to any significant additional (f) dependence. 
It only modifies the 9 distribution with a forward-backward asymmetry term of a few percent. 

A cos20 asymmetry analogous to the Collins effect asymmetry arises in e~^e~ — > t^t^ —* TT^DrTr^VT (and 
other such weak decays). Using the known r weak decay parameters this asymmetry can be calculated entirely 
within the electroweak sector of the Standard Model. A comparison of the data to this Standard Model result 
could provide a good way to check for systematic effects. 

Based on all these considerations it appears feasible to arrive at an actual measurement of the Collins effect 
and an extraction of the Collins fragmentation function to a reasonable extent. Based on arguments in favor of 
the universality of the Collins fragmentation function for the processes e"'"e~ —f hih2X and SIDIS, a simultaneous 
fit to the Collins effect asymmetry data can be and already has been performed, leading to a first extraction of 
transversity [8| . This demonstrates clearly and explicitly the merit of the Collins effect asymmetry measurement 
at BELLE. 
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